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AtomCT 


Heat  tran>.:ftfr  in  boundary  layer  flow  at  high  relocities  and  high 
tenoerstores  is  discussed,  and  engineering  relations  for  its  calculation 
are  presented.  A  similar  discussion  is  the  subject  of  I«aDC  Technical 
Fleport  5I4-7U,  published  in  19SU*  Aeronautical  end  space  flight 
applications  have,  in  the  meantine,  pushed  the  conditions  under  which 
heat  transfer  knowledge  is  iiopcrtant  to  much  higher  velocities  and 
temperatures.  This  need  for  information  created  an  intensive  research 
effort,  wiiicn  is  reflected  in  the  large  number  of  publications 
contained  in  appendix  III  of  this  report.  h(f  this  research,  the 
influence  of  very  large  temperature  aif Terences,  of  dissociation  and 
ionization,  of  low  densities,  and  of  chemical  reactions  on  heat  transfer 
has  been  well  established  in  its  general  features.  In  addition,  new 
cooling  HMthods,  like  transpiration,  ablation,  or  film  coding  bave 
been  studied.  The  reference  enthalpy  method,  introduced  in  dAhv  IH 
5ii-70,  and  well  established  in  the  meantim,  is  extended  end  generalised 
in  this  rorort  to  include  the  additional  effects  nentioned  above.  This 
leads  to  sinple  relations  by  >du.ch  heat  transfer  can  be  calculated  with 
an  accuracy  which  should  be  sufficient  for  engineering  purposes.  Ihese 
relations  additionally  hare  the  advantage  that  most  of  then  are 
applicable  to  any  gas  or  gas  mixture;  therefore,  not  only  to  the 
reentry  problem  but  also  to  heat  tronrfer  problens  arising  in  rockets 
or  other  propulsion  ^sterns  and  in  other  applications. 
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1.  INTBc^OiJGTIUlj 


Klectronic  computers  have  radically  influenced  the  extent  and  Iftnd 
of  arelyticel  work  which  is  today  pe;’for»ed  in  connection  with 
engineering  deTelopnsnts.  inalysea  which  a  few  years  ago  were  st;!!! 
considered  much  too  tedious  azid  rime-consuming  can  be  done  today  in  a 
comparatively  short  time.  However,  regardless  of  the  tremendous  heJ^: 
which  the  computers  offer  to  engineering  designs,  there  still  exists  a 
definite  tned  for  simple  relations  which  can  be  used  in  hand  calculations, 
ouch  x^ations  auwr  needed  for  preliminary  design  work,  or  they  are  also 
useful  for  more  extended  calculations  on  electronic  computers  if  ths 
heat  trani-fer  is  onl^  a  small  part  of  the  total  program. 

The  purpose  of  the  present  report  is  to  present  a  survey  of  the 
field  of  heat  transfer  from  gases  to  surfaces  at  high  velocities  aixl 
high  temperatures,  situations  as  they  occur  in  aeronautical  enginoerli^ 
in  the  development  of  high  speed  vehicles  like  missiles  and  satellites 
or  in  space  flight.  dii:p}.e  relationstdps  are  presented  or  developed 
with  which  heat  tranrfer  can  be  calculated  with  an  accuracy  which 
necessarily  is  limited,  but  which  should  be  sufficient  for  normal 
engineering  design  purposes.  The  discussion  of  the  field  and  the 
relations  which  will  be  presented  are  kept  general  sc  that  they  can  be 
applied  not  only  t;.  external  heat  transfer  on  vehicles  moving  through 
the  atmosphere  of  oui*  earth,  but  also  to  other  gases  at  other  temperature 
and  pressure  conditions.  In  this  way  they  can,  for  instance,  be  used  to 
make  calculations  on  beat  transfer  in  atmospheres  as  tney  exist  on 
otaer  planets  or  to  heat  transfer  in  the  ix>ssles  of  rockets  through 
which  combustion  gases  of  various  comoosition  nre  exhausted. 
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A  similir  sunrey  '^ua  ;;iid«  iu  1^34  and  has  been  publishoU  under 
the  title  ":."irvey  on  Heat  Transfer  at  High  ojpeeds”  as  Wkix;  TR  5U-70. 

The  present  report  can  be  considered  as  an  extension  of  this  one 
considering;  the  chantes  in  the  c-  nditiont  as  they  have  occurred  in  new 
engineering  development  and  of  advances  in  our  knowledge  reported  in 
the  int?r.-ening  years. 

tiesearch  in  the  field  of  heat  transfer  as  it  is  contained  in  this 
renort  has  been  very  intensive  during  recant  years.  The  Appendix  to 
this  report  contains  a  collection  of  pertinent  literature,  and  a  special 
atten^'t  v;us  node  to  make  thir  listi^ig  complete  and  up-to-date. 


2.  aiAhof  BAUhcas 

special  cooling  methods  like  mass  transfer  cooling^  transpiration 
coding,  or  ablation  cooling  are  frequently  used  in  raeent  engineerir^ 
applications.  Parameters  which  'ieseribe  heat  transfer  to  surfaces  cooled 
by  «;uch  methods,  like  recoveiy  factors  and  heat  transfer  eoeffieients, 
ore  defined  in  various  ways  in  the  literature.  Zt  is  therefo’*e  necessary 
to  introduce  in  this  chapter  the  purenetere  which  will  be  used  in  this 
paper.  The  discussion  will  start  with  an  enez^  balance  on  a  solid 
surface,  and  then  proceed  to  surfaces  cooled  by  transpiration,  film, 
or  ablation  cooling. 

Figure  la  indicates  a  solid  wall  with  a  surface  v.  An  energy 
balance  can  be  set  up  for  this  surface  expressing  the  fact  that  the 
amount  of  heat  leaving  the  surface  per  xmit  tine  on  t.he  fluid  side  must 
be  equal  to  the  amount  of  heat  entering  the  surface  on  the  solid  side. 
Heat  will  leave  on  the  fluid  side  by  conduction  and  by  radiation.  Heat 
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soli::  side  b/  conduotir-n. 


-r-.T-er  or. 
l.bis  I’aci  reads: 

-  k 


/vTi  equaLi  n  which  states 


'.4^ 

J  V 


Q  y 


{2-1) 


a  « 


In  this  eiuation,  k  Iniicates  the  heat  conductivity  on  the  gas  side, 
k  the  heat  coiiductlvity  of  the  solid  material,  T  is  the  gas  temcerataire, 


r  the 


temperature  of  the  solid  wall  material,  indicates  heat 


transferred  from  the  well  to  the  sviTounding  by  radiation,  and  q^  the 
boat  conducted  from  the  interior  of  the  solid  wall  tc  the  surface. 


y  i.ndicate:  the  direction  normal  to  the  wall  surface  on  the  fluid  side 
am  on  the  solid  side.  Beth  are  counted  positlYe  in  the  direction 
away  from  the  surface. 

The  temperature  which  the  wall  surface  assumes  when  ths  hsst 
conduction  into  ths  interior  and  the  hsat  loss  by  radiation  are  both 
aero,  is  called  its  recovery  temperature,  T^: 

ly  *  Tj.1  iy  •  iy»  when  On  5  0,  q^^  s  0  (2-11) 

1  indicates  the  enthalpy  of  ths  gas.  ‘ilM  hsat  flux  by  conduction  In 
the  gas  to  the  wall  surface  is  conventionsUy  expressed  through  a  heat 
transfer  coefficient  h,  wnioh  is  defined  by  the  following  sq[uationt 


In  Reference  11.0  it  has  been  demonstrated  that  it  is  sdvsntageous  to 
re-deflm  the  best  transfer  coefficient  by  referring  it  to  snthslpiss  1 
instesd  of  to  temperatures.  This  hsat  transfer  coefficient  is  denoted 
by  h^.  In  this  wsy,  simple  relations  whion  navu  been  obtained  for  a 
constant  property  fluid  can  be  v%3ed  to  describe  heat  transfer  with 
good  accuracy  if  the  properties  of  the  fluid,  including  its  specific 
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iwat,  var;-  widely.  a<iua:.ion  (2-3/  indicates  the  .finition  of  the  heat 
transfer  coefficient  referred  to  entnalpiea  on  tne  .-ijiht  hand  side.  The 


energy  balance  as  contained  in  equation  (2-1)  nay  be  re-written  in 
entha]  :;/.es  in  the  following  way: 


hi(i^  -  ij  +  ‘Vad  " 

Mt>  vill  rw)w  consider  the  transpiraticn  cooling  process.  Figure  lb 
indicates  a  porous  wall  through  which  a  coolant  gas  is  flowing.  A  heat 
balance  will  again  be  nade  statir«  the  fact  that  the  sm  of  all  beat 
fluxes  leaving  the  surface  w  of  the  transpiration  coolec  wall  mist  be 
equal  tc  the  sun  of  tr  e  heat  fluxes  traTelling  In  the  wall  towani  its 
surface.  Ihe  following  equation  describes  this  heat  balance: 


w  w 


cw 


(2-5) 


Two  additional  terns  appear  in  this  equation  because  a  stress  m  of 
coolant  gas  continuously  passes  the  wall  surface  w.  Heat  is  carried 
in  the  wall  by  convection  toward  the  surface,  and  on  the  other  side 
heat  leaves  the  surface  by  convection,  ihe  two  correaponding  terns 
are  the  last  ones  on  both  sides  of  equation  (2-5).  A  is  the 


velocity  of  the  coolant  gas,  Cp  is  its  specific  heat,  and  is  the 
temperature  with  which  the  coolant  gas  passes  the  aurfaee  w.  It  My 
be  observed  that  the  coxrrectlve  torws  cancel  on  both  sides  of  the 
equation,  and  that  in  this  w^y  the  equation  siuplifies  to: 


hi  -  V  ♦  Rrad  =  ‘W  ^^-6) 

In  this  equation,  a  heat  transfer  coefficient  has  bean  introducad  to 
describe  the  heat  flow  in  the  gas  to  the  wall  by  conduction.  It  is 
important  to  stress  that  this  definition  of  tbs  heat  transfer  coefficient 
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will  be  used  in  the  present  report.  In  the  literature,  the  heat  transfer 
coefficient  is  sometines  defined  to  describe  not  only  the  heat  flow  by 
conduction,  but  also  the  convective  contribution  c:c  wording  to  the 
foUoutng  equation: 

hi  (1^  -  V  =  -  k  S  ^cw 

It  is  fouivi,  however,  that  the  definition  as  used  in  equation  (2-6)  is 
advantageous  because  the  heat  transfer  coefficient  defined  in  this  way 
depends  on  a  snaller  nuaber  of  paraaKtsrs.  fhe  heat  flux  by  eonduetion 
can  again  be  referred  to  an  enthalpy  difference  Instead  of  a  temperature 
difference.  The  definition  In  equation  (2-6)  will,  therefore,  be  used 
in  this  report. 

For  steady  state  conditions,  it  is  often  advantageous  to  make  a  heat 
balance  net  for  the  surface  of  the  wall  but  for  a  control  volume  as 
indicated  on  the  right  hand  side  of  Figure  lb.  Also  indicated  in  this 
sketch  are  the  various  heat  fluxes  which  leave  or  enter  the  conirol 
volume.  Cne  surface  of  this  volume  nay  coincide  with  the  wall  surface 
w.  Un  the  coolant  entry  side,  the  surface  is  arranged  so  that  it 
includes  Uw  boundary  layer  which  also  exists  on  this  side  of  the  wall 
surface.  In  this  way,  no  heat  transport  by  conduction  will  occur  through 
this  control  surface,  since  it  is  located  cutside  the  boundaiy  layer. 

The  following  equation  states  the  heat  balance: 

''l“.  -  V  ^  W*  “c.  +  V  =  “c 

^cw  ^  equation  indicates  the  enthalpy  of  the  coolant  gas  as  it 
passes  the  wall  surface  w,  and  indicates  the  enthalpy  of  the  coolant 
gas  on  the  entry  side  of  tlw  wall  and  outside  of  the  bomwlarv  layer, 
q^  indicates  a  beat  flow  which  raay  leave  the  control  volume  in  a 
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direction  parallel  to  the  wall  surface  w,  liquatioL  (2-7.>  has  the 
advantage  that  the  various  terms  appearing  in  it  are  relatively  easily 
accessible  measu.'enent. 

In  the  ablation  cooling  process,  the  heal  flux  from  the  hot  gas  to 
t)ie  wall  surface  is  such  that  the  wall  material  sublimes  and  is  gradually 
carried  av;ay  with  the  gas  stream,  ibe  surface  of  tlie  wall  in  this  way 
gradually  recedes  in  direction  towards  the  wall  interior.  Me 
consider  a  heat  balance  for  the  interface  between  the  solid  wall  and 
the  gas— this  means  for  the  surface  which  gradually  recedes  into  the 
solia  material t  Accordirgly  a  mass  flow  occurs  through  this  surface 
whose  magnitude  per  unit  time  and  area  n^y  be  A.  If  the  enthalpy  of 
the  subliming  wall  malarial  in  the  gaseous  state  and  at  the  temperature 
of  the  interface  is  indicated  by  ig>  and  at  tha  same  temperature  but 
in  the  solid  state  by  then  the  heat  balance  for  the  heat  fluxes 
through  such  a  surface  can  be  expressed  by  the  equation: 

-  V  ("57  )  +  (2-8) 

s  ^ 

'Ihe  enthalpies  and  i^  used  in  the  definition  of  the  heat  transfer 
coefficient  h^  will  be  discussed  in  more  detail  later  on  in  the  paper. 

Tne  equatiob  can  be  simplified  to  the  following  fom  if  the  heat  of 

sublimation  i  s  i  -  i  is  introduced; 

JBB  g  8 

-  V  S-ad"^  I 

In  some  cases,  ti:s  solid  wall  material  is  not  sublimed  but  first  melts 
and  then  evaporstes  into  the  gas  stream,  equation  (2-9)  is  still  valid 
for  this  situation  if  i^^  is  interpreted  as  the  heat  of  evaporation. 
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and  when  the  equation  is  used  to  describe  the  heat  bala  xe  on  the  inter¬ 
face  between  the  liquid  film  which  covers  the  solid  wall  and  the  gas. 
ibe  ablation  process  majr  als/^  be  affected  bgr  sone  chemical  reaction 
between  the  gas  and  the  solid  Etaterlal.  In  this  ease,  the  heat  released 
by  the  chemical  reaeticn  has  to  be  introduced  for  This  term  must 

also  include  surface  dissociation  or  surface  reconblnation  if  it  should 
occur  in  the  actual  process. 

A  quasi-steady  condition  often  occurs  in  an  ablation  process  in 
which  the  temperature  of  the  interphase  is  constant  and  a  temperature 
field  which  docs  n~t  char,£e  its  ‘«hape  moves  with  coaetant  velocity  into 
the  solid  materia!.  For  such  a  condition,  it  ia  again  advantageous  to 
write  an  energy  balance  for  a  control  volume  as  indicated  on  the  right 
hand  side  of  Figure  Ic.  One  surface  of  this  control  volume  coincides 
with  the  interface  w  of  the  ablating  material}  the  other  surface  is 
parallel  to  this  interface  but  at  a  stiff icient  distance  so  ttu  i  the 
temperature  increase  occurring  in  the  wall  by  conduction  has  not  yet 
reached  this  point.  The  individual  heat  fluxes  leaving  this  control 
volume  are  again  indicated  in  the  figure,  and  it  is  assumed  that  no 
heat  will  leave  the  volume  in  a  direction  parallel  to  the  wall  surface. 
The  following  equation  then  holds  for  this  heat  balance: 

<»rad  *^g  =  *^so  ^2-10) 

i^  is  the  entlialpy  of  the  solid  wall  material  at  tenpermture  T^. 

If  the  ablation  process  occurs  as  a  melting  and  evaporation,  and  if 
/3  Indicates  the  fraction  of  tho  ablated  material  which  leaves  the 
control  volums  in  liquid  form,  then  the  equation  deseribii^;  the  heat 
balance  is: 
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h(i  -  i)  +  q  j  +  0  m(i  -  i  ;  +  -  f}  y  m  (i_  -  i  }  »  0  (2-11) 

1  w  r  rad  '  g  so  '  f  so 

ic  the  eiiths'oy  of  the  liquid  naterlal  leaving  the  control  volume, 
ihe  balance  equations  contai;.ci  in  this  iection  have  to  be  used 
in  cc>nnec<.ion  trith  other  relations,  for  instance,  one  vhich  describes  the 
heat  corehiction  process  in  the  solid  material  for  the  various  engineering 
calculations,  A  prerequisite  of  their  use  is  a  knowledge  of  the  heat 
transfer  coefficient  ana  of  the  recoverj'  tenperature  or  enthalpy. 

Iheir  discussion  is  the  subject  of  tr.e  following  chapters  of  this  paper. 


3,  Haal  'itiaHSeaB.  tO  Sv-tli)  3Ud?AClti 
31,  SlliaLE  COrlfOIi^tT  GaS 

Heat  transfer  to  a  gas  of  high  tempera tui-e  is  influenced  by  the 
fact  that  the  gas  aa^  be  partially  or  completely  dissociated  or  ionised, 
and  that  the  degree  of  dissociation  changes  throughout  b'.'unuaxy 
layer  when  the  temperature  of  the  surface  eapesad  to  the  flow  is  different 
from  the  tenperature  of  the  gas.  Such  effects  may  also  occur  if  a  gas 
with  a  relatively  low  teaoerature  flows  with  very  high  velocity  over 
the  surface,  because  the  gas  is  then  heated  within  the  boundary  layer 
by  internal  friction,  and  the  increased  tenperature  may  again  causa 
dissociation.  Ihe  influence  of  such  factors  cn  heat  tratsfar  will  be 
discussed  in  the  next  section  of  tiiis  report,  m  the  present  section, 
sucit  effects  will  be  neglected.  It  is  therefore  assuaed  that  we  deal 
with  a  si:igle  conponent  gas,  and  that  no  dissociation  or  reconbination 
in  the  gas  occurs  as  a  consequanca  of  locally  vazying  teeiperatara. 
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31x»  1-lat  -late 


Heat  transfer  to  the  -urface  of  an  object  exposed  to  a  high 
velocity  flow  de  ends  on  the  shape  of  tl\is  object  as  well  as  on  the  flow 
field  to  whi';i  it  is  exposed,  une  geoitetry,  which  has  been  investigated 
very  extensively  because  it  can  be  treated  theoretically  in  a  siaple 
aanner,  is  a  situation  conventionally  referred  to  as  "flat  plate."  kb 
define  as  a  "flat  plate''  situation  one  under  which  a  plana  surface  Is 
exposed  to  a  flow  field  with  locally  unifom  velocity  and  flow  direction. 
It  is  asstuned  that  the  flow  velocity  ancL,  as  a  eonaequence,  the  pressure 
arc  also  constant  along  the  outer  edge  of  the  boundary  layer  which  exists 
along  the  plane  suriace.  It  is  fortlter  specified  that  the  boundary 
layer  is  so  thin  that  the  p  -essure  is  transAitted  without  change  to 
'the  surface  of  the  plate,  ihis  flow  condition  can  be  set  up  experlMa- 
tally  with  good  approxissatlon  in  low  velocity  flow  when  the  Beynolda 
nuii^r  ia  aufflciently  high.  In  supersonic  and  lypersonie  flow,  the 
boundary  layers  axe  usually  thicker  and,  aa  a  consequence,  a  saock 
wave  is  generated  by  the  boundary  la/er  ahead  of  the  plate.  As  a 
consequence,  the  velocity  and  pressure  vary  along  the  outer  edge  of  the 
boundary  layer  and  also  along  the  plate  surface.  Ihe  temperature, 
which  is  also  constant  along  the  outer  edge  of  the  boundary  layer  for 
"flat  plate"  flow,  varies  locally  where  a  shock  is  created  Iqr  the 
boundary  layer.  In  the  present  discussion  these  effects  will  be 
disregarded,  and  the  condition  which  has  been  aentioned  above  of  a 
locally  constant  pressure  and  of  a  constant  tenperature  or  enthalpy 
along  the  outer  edge  of  the  boundary  layer  will  be  specified.  In  this 
way  a  standard  case  is  created  for  idtich  siraple  relationships  describe 
the  heat  transfer.  These  relations  can  also  be  used  to  approxiaats 
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Jfi. 


coac:it,ions  i:;  sj.eno.er  objects  in  supersonic  flow  when  the  region  near 
tne  leading  edge  is  exclucen,  and  wtien  the  actual  velocity,  pressure, 
and  temperature  at  the  uter  edge  of  the  boundary  layer  are  in-oroducod. 
fhenc  may  be  considerably  different  from  the  values  in  the  free  stream. 


3111«  Constant  aurface  lenoerature.  iihren  under  the  situation 
which  has  just  been  defined,  heat  transfer  is  not  uniquely  described, 
it  still  depends  on  the  way  in  which  the  temperature  cn  the  plate 
surface  varies  locally.  In  this  naragraoh  it  will  be  assumed  that  the 
stirface  temi'erature  is  l''cally  constant,  fhe  influence  of  a  local 
variation  of  this  tempera  lAire  will  be  discussed  in  tlie  f-'llowing 
paragraph. 


if  the  fluid  to  which  the  plate  surface  is  exposed  hajt  constant 
properties  (independent  of  temperature  and  pressure),  then  heat  transfer 
as  well  as  friction  are  described  by  swiple  relations  which  were  derived 
for  a  laminar  boumary  layer  a  long  ti  »  ago,  and  whici  have  been 


experistentally  verified  to  a  high  degree  in  the  meantime,  similar 
relations  from  experimental  results  have  been  developed  for  a  turbulent 
boundary  layer,  xbese  relations  are  customarily  expressed  in  dimension- 
less  parameters  defined  in  the  following  way.  The  recovery  temperature 
in  equation  (2-2)  is  expressed  by  a  dimensionless  parameter  r  called 
temperature  recovery  factor  and  defined  the  following  way: 


(3-1) 


Ug  is  the  velocity  aiil  the  static  temperature  at  the  outer  boundary 
layer  edge.  I'he  heat  transfer'  coefficient  ir:  equation  (2-3/  is  expressed 
by  tire  dimensionless  iJusselt  number 
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Nu= 


(3-^) 

where  x  is  t.hs  distarice  from  tne  leading  edge  measured  aluug  the  pxo^e 


surface,  or  altemateljr,  by  another  one  called  ^tanton  number 
St  •  - - - 


(3-3) 


p  c  u 
T  p  e 

^dtere  ^  is  the  density  and  the  specific  heat  at  constant  pressure 
of  the  fluid.  In  some  situations  the  shear  which  the  flow  exerts 
at  the  plate  surface  is  of  interest  as  hoU  as  the  heat  transfer.  Ihs 
local  shear  is  deterainad  hgr  a  dimensionless  parameter  c^  called 
friction  factor: 


"f  • 


tw 


(3^) 


i'or  laminar  flow  conditions,  the  following  relations  express  the 
tcD^rature  recorexy  factor,  the  lusselt  isMber,  Stanton  number,  and 
friction  factor: 


r  *  (3-'^) 

f3r  is  the  frandtl  nnter  of  the  fluid. 

1/3  j— 

Ku  •  0.332  (ftp)  (3-6) 

n  tux 

Be  <■  3  is  the  Reynolds  number  and  M  tbs  dyztamie  viscosity  of 

r 

the  fluid. 

2/} 

St  =  0.332/{ftp)  rS  (3-7) 

s  0.66)4/  (3-8) 

It  can  be  ssea  that  tbs  following  rel.-.tion  exists  between  the  Stanton 
aa^or  and  the  frietics  factor: 

-  2/3 

St  s  (  cj/2)(Pr)  (3-9) 

For  turbulent  flow,  the  following  relation  was  found  from  experiments 
to  dsaeribs  the  temperature  recoveiy  factor: 
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V 


(3-10) 


r= 

The  locfl  frictioii  factor  has  been  described  by  Yarious  iaaJytic 
»^x:<rcssic  .s.  "  arnple  one  is,  for  :.nstance.  the  equation  developed  by 


we.'ju.!  ‘.i-JiTinow; 


(J-Ii) 


It  has  also  been  foual  that  the  relation  between  dtanton  nuatber  and 
friction  factor  in  turbulent  flow  is  practical../  the  sane  as  for  laminar 
flow  conditions,  sc  that  equation  13-9)  hold.s  for  turbulont  flow  edso. 
i'ron  this  relation  H  follows  that  the  dtanton  number  for  a  turbulent 
boundary  layer  on  a  flat  p.Iate  Is  described  by  the  following  equation: 


(3-12) 


■or  real  c'-ses,  it  has  been  fosind  that  the  properties  are  not 
consent  but  vary  with  ta.j.>rature.  ihc;.  can,  Imwever,  ■'e  considered 
as  independent  of  pressure  os  long  as  excessively  high  pressures  and 
dlssc  slat  ion  are  er.cluded.  Under  this  condition,  and  as  long  as  the 
specific  }*eit  Cp  can  be  coi:siticred  constwit.,  it  has  been  established 
in  reference  110  that  the  relations  jlven  In  the  preceeding  paragraphs 
fo:  u  constant  property  fluid  still  describe  the  actual  heat  transfer 
ar/^  f-ic-lon  witn  good  accur.-.cy  when  :tl»e  properties  xpoearlng  in  all  of 
these  relationships  arc  intr-duced  into  trio.S2  relations  at  »  reference 


tenparaturt  described  by  the  following  equation; 

T*  s  Tg  0.!;(i^  -  TgJ  +  0,-?2(Tj.  -  T„) 


(3-13) 


Tho  secTd  ten  v.i  t-hic  equation  can  be  expressed  by  i-jach  numtar 
Ma^  exJ.sting  in  the  flt-rf  a.  the  outer  edge  of  the  boundary  layer.  The 
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equa'-.i':;!  then  '.ou.ca  the  fcliowing  fcm  which  is  sonetr-.js  more  conven¬ 


ient-  for  an  evaluation: 

f*  t  'i’g  +  O.JCf,,  -  Tg)  0.22r  (3-lli) 

^  is  the  ratio  of  specific  heats  at  constant  pressure  and  voltine. 

Vftien  the  tenperature  variation  within  the  boundary  layer  is  so 
lai-ge  that  the  specific  heat  of  the  varies  considerably,  then  the 
better  procedure  is  to  base  heat  transfer  calculations  on  enthalpies 
instead  of  on  temperatures,  in  txiis  case  tne  definitions,  as  given  on 
the  xdght  hand  side  of  equations  (2-2)  and  (2-3),  are  used.  The  heat 
transfer  coefficient  hj^  is  •3:q>ie'oeJ  in  dimensionless  form  as  a  Itasselt 
number  or  a  Stanton  ooiber: 


NUi 


»  Stj^ 


.  *ii 


(3-15) 


Ic  " 

The  relations  (3-5)  to  (3-12)  can  again  be  used  to  describe  these  nee 
parameters.  Ibe  properties  appearing  in  these  equations  are  now 
introduced  at  a  reference  enthalpy  instead  of  a  reference  temp'^rature. 
The  following  equations  describe  the  reference  enthalpy: 

i*  s  o.5(i^  -  ig)  +  0.22(ij.  -  ig)  (3-16) 

«  -  1  o 

i  s  ig  +  0.5  (Iv  -  ie>  +  0.22r  — (Kag)‘‘  i,  (3-17) 


fna  above  procedure  based  on  a  reference  enthalpy  and  on  constant  proper^ 
relations  is  tudv  in  widespread  use  and  is  usually  referred  to  as  reference 
enthalpy  method.  Its  accuracy  has  been  checked  by  oo«q>arison  with  the 
results  obtained  fay  bointdaiy  layer  solutions  for  a  laminar  boundary 
liyrer  and  for  air  as  the  flowing  mrdlua.  Agreement  within  plus/ainus 

15  has  been  found  for  stream  temperatures  between  iiX)  and  300°R,  for 
vail  temperatures  between  kjO  and  2000°R,  and  for  Mach  numbers  up  to 

16  (Ref.  110).  It  will  be  seen  in  a  later  section  that  the  method  still 
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is  valid  even  i:  dissociation  or  chemical  reactions  occur  within  the 
boundary  layer.  For  ♦his  r'oas-n,  the  discussion  In  this  report  will 
ma3pj.y  be  based  on  the  reference  enthalpy  method.  Agreement  is  also 
good  with  measunsments  on  turbulent  heat  transfer  in  the  l^ach  rujoher 
range  which  have  been  covered  up-to-date.  It  is  also  reasonable  to 
expect  that  the  relations  hold  witli  good  accuracy  for  a  gas  diffe?.'er.t 
from  air. 

3112.  i/ariable  Surface  Tenoerature.  ihe  various  methods  ^lich 
have  been  developed  to  calc’ilate  heat  transfer  to  a.  surface  with 
locally  varying  temperature  have  -Iready  been  discussed  in  reference 
liO.  A  simplified  procedure  which  allows  the  calculation  of  heat 
transfer  to  a  wall  with  prescribed  varying  temperature  has  been  developed 
in  the  meantine  (Ref.  113).  '^'he  procedure  approximates  the  actual 
temperature  variation  by  a  succession  of  straight  lines  snd,  in  this 
way,  simplifies  the  calculation  to  a  sumation  process*  It  will  bo 
discussed  in  connection  with  Fig.  2.  This  figure  indicates  an 
arbitrary  variation  of  tno  temperature  potential  s  wldoh 

is  detenaining  the  heat  transfer  according  to  equation  (2»3J  •  i’he 
surface  of  length  x  is  subdivided  into  a  number  of  partial  lengths 
dx  of  equal  dimensions.  The  temperature  potential  at  the  positions 

*o»  *2»  ^3***^n  ^  following 

equation  con  be  used  to  osloulate  the  looal  heat  flux  qu,oo 

time  and  area  at  the  position  x^^,  and  a  very  similar  equation  allows 

I 

calculation  of  the  total  heat  transferred  along  the  plate  surface 

between  location  U  and  xi 
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'1«,C0  *  *’i50  1^0  - 

'd'g  -  2  iJ'^  ~  2 13;  ..V  -  2^^  _  ^]}  (>.16) 

Vco=^i»o^{’^o''^‘^^n-  -B'(4Vx)  i(2n-V^^- 

,^o-2>,  -21^2--..  (3-19) 

h.  is  the  local  beat  transfer  coefficient  and  h.  the  aTerage  me 
iso  i*o 

for  the  saiiifa  flow  condition  but  for  a  constant  wall,  tenperature,  b  is 
the  width  of  the  plate  noznal  to  flow  direction,  xhe  following  equation 
can  be  used  to  calculate  the  isotbenoal  heat  transfer  coefficient: 

»i»  = 

The  exponent  n  is  equal  ),$  for  laminar  flow  and  equal  0.8  for  turbulent 
flow.  The  constants,  a,  B,  C,  a',  B',  0*  contained  in  these  equations 
are  listed  in  Table  1«  The  last  en/jaticn  repres?-nts  f'C  average 
heat  transfer  coefficient  when  is  replaced  by  C*» 


TABLE  1 


IB 

A 

B 

c 

A* 

B' 

C 

l>aninar  0 

d.o95 

0,o90 

0.332 

0.696 

O.U32 

0.661* 

vu 

0.837 

0.635 

o.inz 

1/2 

0.81P 

0.572 

0.1*69 

1 

0.792 

0.53? 

0.560 

Turbulent  0 

0.951 

0.117 

0.0296 

0.962 

O.U78 

0.0370 

The  condition  n  *  0  aptiO^ies  to  flat  j^ate  flow.  The  values  for  a 
paranctars  different  from  0  aptly  to  flow  with  nressure  gradients  and 
will  be  discussed  later  on  in  this  report. 

Scisting  methods  for  the  calculation  of  heat  transfer  with  a 
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Lamtnar 

turbuhift 


Ratio  of  Nusselt  number  cn  surface  with  varyinp  temperature 
to  Vusselt  number  on  isothermal  wail,  (fron  Ref.  113) 


locally  varying  surface  tenpcrature  assune  a  fluid  with  constr.nt 
properties  or  a  fluid  for  which  the  oroduct,  density  tiaes  viscosity, 
is  a  ccnstajit-  i.he  presented  ap.nDxiiitation,  therefore,  is  based  on  the 
sane  ass'imption.  It  checks  the  results  of  ca..culation5  oerfomed  with 
the  other  aethods  within  a  few  per  cent. 

For  a  fluid  with  variable  properties,  it  is  prt^sed  to  coxntine 
the  nietliod  presented  here  with  the  reference  enthalpy  riethod  of  the 
preceding  section  by  introducing  properties  at  the  reference  enthalpy 
into  all  of  the  relations  api'carir^  in  the  equations,  rcqperiiaents 
under  conditi.  ns  where  the  wall  temperature  v..ries  widely,  aiw  where 
the  temperature  differences  in  the  boundary  layer  are  sufficiently 
large  are  r.oL  extersiTe  enough  to  check  tne  accuracy  of  a  calCttlatioB 
made  with  the  proposal  just  .uintioned. 

Ths  influence  of  a  wall  tenpr-rature  variation  on  heat  transfer 
becomes  significant  only  v.-hen  the  temperature  variation  is  large  relative 
to  the  te*i:'«rature  i>ctential.  a  first  estimate  of  the  itulueacc  of  a 
wall  temperature  variation  and  a  decision  whether  it  has  to  be  considered 
in  a  dete:riination  pf  t'ue  heat  transfer  can  be  made  ..’ith  the  help  of 
fig.  3«  In  -.his  figure,  t;»  actual  .iusselt  number  dirtded  by  the  iiaseelt 
raimber  at  an  isothemal  surface  witii  ti:e  sane  local  temperature 
difference  is  plot.ed  over  a  oarameter  ^  .  ihe  wail  temperature 
variation  is  assumed  to  follow  an  exponential  laiw  as  indicated  in  the 
figure  idtri  the  value  ^  as  exponent.  Vlie  curves  with  the  parameters 
m  ■  0  apply  to  flat  plate  conditions  and  curves  with  a  finite  m  value 
tr  flows  with  pressure  gradients  wiiich  will  be  discuS'-ad  later  on.  It 
may  be  obsen'ed  that  the  influence  of  a  tent-erature  variation  is 
considerably  lfi.*ger  for  la-iinar  f).ow  conditions  than  for  turbulent  flcv. 
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ror  turbulent  flow  it  is  only  nacessarj'  in  very  rare  cases  to  make  a 
correction  for  tiie  temperature  variaticn  cr  : heat  ti'  -.osfer 
coefficients. 


112,  Blunt  Ctolects 

aarly  develo-vwnts  in  high  speed  fli^.t  were  concerned  witli  objects 
of  very  slender  ihape.  fhis  was  done  in  ordei  u)  decrease  the  drag  and 
the  power  required  to  aove  such  objects  tnrough  the  ataospbere.  In 
new  uevclcpncnts  and  in  astronautics*  on  the  other  hand,  interest 
usually  centers  arouni  blunt  objects,  becattse  dreg  is  of  minor  inporMUice. 
donetines  drag  is  even  desired  for  a  reduction  of  aerodyiramic  heating, 
duch  an  object  flying  with  h/pci^s^nic  veloeitieB  finds  'sonditions  in  its 
surroundings  as  sketched  in  i-lg.  a  boiudaiy  layer  envelopes  the 
object  as  a  thin  sheet.  Outside  of  the  boundary  layer,  e  shock  is 
created.  In  the  bae<<  of  the  object  the  flow  usually  separates,  so  that 
the  region  behind  the  r'earward  surface  is  filJed  with  a  :low  containing 
strong  irregular  vortices,  i'be  density  in  this  region  is  usually 
quite  low. 

In  moving  through  the  shock  wave,  the  air  loses  most  of  its 
kinetic  energy  ant;  converts  it  into  internal  energy.  As  a  consequence, 
the  iiyer  oeU-eer.  the  shock  and  the  outer  edge  of  tJie  boundary  layer 
is  very  hot.  fhe  conditions  at  the  outer  edge  of  the  boundary  layer 
which  are  required  for  a  calculation  of  heat  transfer  to  ti^ie  surface 
can  be  obtained  with  good  accttracy  from  the  80,>called  Newtonian  flow 
approxinaiion  for  nach  nusters  treater  than  approxinstely  !*.  It  is 
found  that  in  this  care  the  shock  wave  is  quite  close  to  the  surface 
of  the  object.  Correspondingly,  ths  velocity  cn"5poner.'-  noraal  to  the 
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■  I'.j  •••  ax_  It.,  ft  .to.ienvu.".  ccn"i;ei:it:on  then 

t,  i  n'  r-.lai.i.-.-i  betveen  t.':e  pressure  Pg  oehind  the 

srecic  vecu.ii  tc-  •-r-o  rr  .-sure  at  tne  -al.  surface  acccrding  to 
be.t.'iJ.:;-..-  layer  trecry;  ani  the  pressure  Pao  and  velocity  u^  in  the 
upstrea.  area  tar.ead  of  tne  sn  cki: 

^  C 

?e  ■  ?«>  =  CO  008*^6 

e  is  the  an,-l<*  bet-.eeii  tne  s^rf^ce  nor.-.-ial  aj>d  the  upstrea-t  flow 
directirn.  ihe  velocity  and  tenperature  alonj;  .e  outer  edge  i.-'  the 
bo  .adary  jayer  c.in  then  be  found  -.-ith  the  assurption  that  the  flow 
expands  outsiac  tne  bounaar;/  layer  isentrotically  frr:i  the  pressure 
at  stagnation  p  int  to  the  local  pressure  p^.  inis  intorccuces  an  error 
which  becuites  larger  witn  increasing  distaiice  fron  tne  stagnation  point, 
bee;  use  t;.e  air  enterii^  the  boundary  layer  at  sodc  distance  from  the 
stagnation  ooint  actual  y  crossed  the  shock  wave  at  a  location  where 
the  shock  is  inclined  towards  tl«  upstream  flow  direction,  .fus  error 
is,  hciwe-.-er,  usually  snail  for  blunt  objects. 

3121.  ..tacnation  rlow.  In  a  limited  region  close  to  the  stagnation 
point,  it  is  founu  that  the  flow  velocity  increases  li.. early  with 
distance  fron  the  stagnation  point,  ihis  region  is  referred  to  as 
"stagnation  point  region",  and  heat  transfer  in  this  region  can  be 
determined  by  exact  solutions  of  the  boundary  layer  equations,  ihe  flow 
in  such  a  region  is  almost  exclusively  laminar.  The  dimensionless  para¬ 
meters  describing  heat  transfer  in  this  region  are  agnin  described  by 
the  equations  i3-15>.  ihe  tern  x  now  indicates  the  distance  from  the 
stagnation  point  measured  along  the  surface  of  the  object  in  flow 
direction;  u^  is  the  local  velocity  of  the  fluid  at  location  x  and 
just  outside  the  boundary  layer.  The  recovery  enthalpy  5  .  in  the 
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stagnation  point  r«gion  is  prsfticaliy  oqual  to  ths  tot«l  gas  snthslcj 
in  toe  upstrean.  in  stagnation  point  flov,  ons  has  to  distinguish 
betwc-^  two*dinensioual  flow  as  it  occurs,  for  instance,  around  tt» 
front  part  of  a  wii^  and  bstMoen  rotationally  tqrametric  flov  as  it 
usuaUjr  exists  around  the  nose  of  a  nisslle.  the  first  situation  will 
be  referred  to  as  "plane  stagnation  flow"  and  the  second  as  "rotationalljr 
synaetric  stagnation  flow",  tor  a  constant  property  fluid  and  a 
lanirar  boundary  layer,  the  following  n;lations  describe  heat  transfer 
for  these  flow  coniitioos: 

plane  stagnation  flow  0.570(1^)  *^V^t  (3-21) 

rot.  sfm,  stagnation  flow  0.763(Pr)^*^  ll^c  (3*22) 

3ms  first  one  of  thct^e  relations  was  obtained  by  bquire.  ibe  second 
one  follows  fit>a  the  first  expression  by  ilsngler’s  transfomstlon.  It 
was  also  directly  calculated  fay  N.  difaulkin.  For  a  fluid  with 
variable  properties,  the  reference  enthalpy  aetbod  was  again  found  to 
describe  heat  transfer  conditions  with  sufficient  engineering  aeevrecy. 
Ihis  was  pointed  out  fay  n,  innig  (Bet'.  296  and  297)  fay  conparisoo  of 
results  obtained  '..-ith  this  aetbod  with  exact  faouadsiy  layer  solutions, 
it  will  be  cnee  acre  deaonstreted  hereby  e  ccapsrtson  of  the  reference 
enthalpy  aetbod  with  a  relation  which  has  been  obtained  by  Fay  and 
Hidden  (Ref.  12U).  The  x'ollowlng  relation  is  contsinsd  in  Bafcrcnce 
12lt  for  a  constant  iYandtl  anaber: 


s  O.b?  ( 


Vw/^  w 


(3-23) 


Tbs  Indsx  w  in  this  equation  indicates  that  properties  arc  to  be 
introduced  at  the  conditions  as  they  exist  at  the  vrll.  aoc.  tdii«  to 

^  According  to  H.  hichbom 
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the  reference  enthalpy  method,  the  properties  have  to  be  introduced  into 
the  constant  pronerty  relations  at  a  reference  condition.  For  either 
plane  or  rotatlonally  syintrie  stagnation  flow  the  following  relation 


I 


exists: 


Mu 


*  C  (FrJ 


Xhe  constant  C  in  this  relation  has  different  values  for  plane  or 
rotatlonally  sjenetrlc  flow,  xhe  reference  enthalpy  *t  which  the 
properties  hare  to  be  introduced  is  given  by  equation  (3-16).  In  the 
■ei^iborhood  of  the  stagnation  point  the  velocity  is  very  Mali,  and 
the  difference  between  the  recovery  enthalpy  i^  and  the  static  esthtO.py 
i^  is  nsgliglble.  Ocarrespondingly,  aquation  (3-l6>  can siopllfied 
to  the  following  expression: 

4I-  s  0.5  a  t  -~  ) 


Xhe  ratio  of  ausselt  to  aquare  root  of  dqmolds  oowber  rvr<.9sc';  at 
refereace  coadi;,ion  can  he  written  in  the  following  way: 


<2*  „•  iTv  FfJ?  TO 

If  both  paraKters  are,  on  the  other  hand,  deternlned  at  wall  conditions. 


then  the  following  condition  holds: 


Tv  »fw<*w  W 


(3-21i) 


a  conparisoa  between  the  equatione  (3*23)  and  (3-2U}  is  aada  in  Fig.  5« 
Ihe  ratio  atraaa  and  at  wall  conditions,  as  used  by  fby  and 

Riddell,  is  shown  in  the  lower  carve  of  this  figure.  3ie  expressions 
presented  by  the  two  upper  curves  are  the  ones  appearing  in  the  eqaatiocs 

Tikvc  re  .'9^2h  25 


V, 


{}~23)  and  (3-i^li^.  'i'he  fact  that  these  two  curves  alrao  t  coincide 
indicates  that  the  reference  enthalpy  method  leads  to  practically  the 
same  result  as  Fay  -  Riddell's  equation. 


3122.  Swept  leading  adges.  A  considerable  rsduction  of  recovery 
temperature  and  heat  transfer  on  the  stagnation  line  of  a  cylindrical 
object  occurs  when  the  cylinder  axis  is  inclined  under  an  angle  smaller 
than  90  degrees  toward  the  upstream  direction,  a  shock  surrounds  the 
front  part  of  such  a  cylinder,  wuen  the  upstream  velocity  is  supersonic. 
At  the  stagnation  line,  the  velocity  component  normal  to  the  cylinder 
surface  is  reduced  t:  zero,  fhe  velocity  component  parallel  to  the 
cylinder  axis,  however,  is  maintained  outside  the  boundary  layer.  Its 
value  may  be  u^.  this  velocity  comp.'nent  is  maintai;«d  in  the 
"stagnation  flow  region",  whereas  the  velocity  component  parallel  to 
the  surface  but  normal  to  the  cylinder  axis  increases  again  linearly 
with  distance  x  (measured  now  in  a  plans  normal  to  the  axis). 
for  such  a  flow  situation  have  been  reported  in  the  literature.  For  a 
constant  property  fluid,  beat  transfer  coefficient  and  recovery  iac-tor 
can  be  again  calculated  from  equations  (3>21)  and  (3-^)  'when  the 
recovery  factor  is  defined  in  the  following  way: 

-  fat 

r  s  — - - 5* 

^'o  *  ^st 

Xq  is  the  total  temperature  in  the  gas  outside  the  boundary  layer, 

and  ^indicates  the  static  temperature  at  the  stagnation  point  outside 

St  2 

the  boundary  layer  (Xq  -  i^e  ).  Xhe  heat  transfer  coefficient 

2c 

is  defined  with  the  temperature  difference  For  a  single 

component  gas  with  variable  properties,  it  is  demonstrated  in  Ref.  itl7 
that  equation  (3-23)  (converted  to  plane  stagnation  flow)  represents 
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tbe  results  of  bojodery  l^yer  solutions  within  in  the  range  of 
parsMters  which  is  of  practical  Interest,  provided  the  heat  ^ronsier 
coefficient  and  recovery  factor  are  based  on  enthalpies  instead  of 
temperatures.  It  is  argued  in  Ref.  Ul?  that  this  should  apply  also 
when  dissociation  and  recoabination  occur.  Kroa  tbe  coaparison  in 
Fig.  it  is  evident  that  th'i  reference  enthalpy  method  &od  sq-ustlone 
(3-5 J  and  (3-21)  describe  beat  transfer  to  the  stagnation  region  of  a 
yawed  oyllnder  with  the  sam  accuracy. 

3123.  Arbitrary  Stapea. 

flow.  Tbe  flow  over  a  tlunt  object  is  connected  vlth  pressure 
gradients  along  the  surface  which  influence  the  develoi«ent  of  the 
bouiilaty  layer  and  of  the  friction  and  iksat  transfer  at  the  surface, 
an  exact  solution  of  the  differential  equations  describing  bou;)dary 
layer  flow  is  under  such  cireunstances  very  tedious,  because,  even  frr 
plane  or  rotationally  sywMtric  flow,  it  mans  the  solution  oi  partial 
differential  equations  in  two  independent  variables.  As  a  cnnwquenes, 

W  * 

nacgr  approximte  procedures  have  been  proposed  which  eircuavent  this 

A 

difficulty.  Tnsse  methods  have  to  be  evaluated  according  to  the  ease 
with  which  they  lead  to  results  and  accordii^  to  their  accuracy.  Una 
which  usually  is  especially  sinole  is  based  on  the  asroaption  of 
'^oeal  siailarity'*.  This  msthod  was  proposed  for  the  first  tlm  by 
Faltaer  anl  ^csn  (V.  K,  Palkner  and  5.  W.  Skan:  rhU.  Hag.  1^  (l93i), 
865*  ieron.  Res.  Conn.  dep.  a.  ‘-wn.  13lit  (3.931))  for  the  calculation 
of  akin  xVlction  and  by  cckert  and  Drewitz  (d.  dckert  and  0.  drewits: 
Luftfahrtforschung,  1^(19lt2),  189)  Tor  the  calculation  of  heat  transfer. 
Haogr  extensions  have  oeen  described  since  that  tine.  Xhs  discussion 

i? 
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in  this  section  will  mainly  follow  a  procedure  sut^ested  by  L.  Lees 
(Ref.  430',  because  it  is  relatiyel;-’  simple  and  can  be  applied  'or  flaiia 
as  well  as  for  rotaticnally  symetric  flow.  The  method  will  be 
somewhat  modified  to  incorporate  the  reference  enthalpy  procedure. 

The  differential  equations  which  describe  steady  laminar  boundary 
layer  flow  of  a  medium  with  variable  properties  are; 

T;  4?  ^3-25) 

^(u 4f + » 4?^  =  - ^  4v  ^  ^3-^^ 

?  4;  “3y^  ^  4?  4?^  +  4?  0-27) 

The  equations  are  written  in  a  coordinate  system  as  indicated  in  Fig.  6* 
The  symbol  s  indicates  the  distance  measured  in  the  flow  direction 'alon( 
the  surface  from  the  stagnation  point,  y  the  distance  from  the  surface, 
and  r^  the  distance  of  the  surface  point  under  considers  lor.  from 
the  axis  of  rotation.  The  index  e  indicates  conditions  at  the  outer 
edge  of  the  boundary  layer.  X  is  the  symbol  for  the  total  enthalpy  in 

u*2 

the  flow  (Xg  *  ig  +  .  The  equations  in  the  above  for»  are 

applicable  for  plane  as  well  as  for  rotationally  sywwtric  flow.  For 
the  first  situation,  n  has  to  be  set  equal  to  zero;  for  rotationally 
symetric  flow,  n  is  equal  1.  It  has  been  shown  by  Oorodnitzyn,  L. 

Howarth,  C.  H.  Illingworth,  and  A,  Stewartson  that  the  equations  can  be 
brought  by  a  traasfomi.r.tion  of  coordinates  into  a  form  which  closely 
resembles  that  for  incompressible  flow  of  a  fluid  with  constant  propertias. 
In  addition,  W.  Mangier  demonstrated  that  the  equations  for  rotationally 
symmetric  flow  can  be  made  to  assume  the  form  of  those  for  a  plane  flow 
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by  another  transformation.  L.  Lees  combined  both  transiormations  by 
the  choice  of  the  following  independent  Toriabless 


2a  ^ 

,  r  da 
e  o 


(3-26) 


2^^  d. 

J  Vo 


if,  in  addition  the  stream  function  which  is  introduced  to  eliadnate 
equation  (3-25)  is  written  in  the  following  fora 

f(->^) 

thesi  the  above  system  of  bounaary  layer  equations  transforms  to: 


(Cf ' ' )  •  +  ff '  •  ‘  P  [  -  (f ')^  r  0 

..  2  ^ 


(3-30) 


(-^  e')'  +  fg'4-~[2C(l  -  -^)f'f"]'aO  (3-31) 


The  following  notations  have  been  used  in  this  eqiiation: 

E  5  -r  »  c  «  1‘A'  ,  (3  s  - 

*  •  V  d» 


(V32) 


*  priao  indica*,ed  differentiation  towards  .  'fbe  assumed  fora  of  the 

stream  function  y  on  irtiich  the  derivation  of  the  equations  (3-30)  and 

(3-31i  were  based  restricts  the  boundary  conditions  to  which  this 

equation  can  be  applied..  A  detailed  discussion  of  this  point  is 

contained  in  the  original  paper  by  Lees  (Ref,  438).  Specifically,  the 

Kiethod  of  local  similarity  postulates  that  the  local  flow  conditions  on 

ary  point  of  the  surface  of  a  body  wi+h  arbitrary  shape  correspond 

u  2 

with  good  aporoximation  to  those  of  a  flow  with  0  and  being 

'  21^ 

constant  and  having  values  corresponding  to  the  actual  local  condition. 
Solutions  of  the  boundary  layer  equations  for  this  condition  are  avail¬ 
able  for  certain  cases.  For  a  vanishing  Mach  number  which  causes  tbs 


•WADC  TR  59-624 


I 


97 


OA 


- 1 

OA 
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A6 


Fyg,7 

"Jtjcnsiciilesr  enthalpy  gradient  on  corled 
surfaoer  (I'rca  213) 


3i 
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last  texm  in  equation  (3-jtii  to  disap  ^ear,  ttoy  hare  been  worked  out  by 
a.  sckert  and  d.  drown  cusi  associates.  For  a  gas  vitii  rrandtl  •  1 
and  =  constant,  solutions  have  been  obtained  by  S,  Levy  (Joum. 
Aeron.  Scl.>  2j^Cl9$b ) .  as  wall  as  by  Cohen  and  Reshntko  (Ref.  66). 

For  tl:  =  1  and  s  constant,  the  aquations  in  rafaresca  66  becona 
identical  to  aquations  (3-30)  and  (3-31)  idtan  the  paraaetar  is 
based  on  Mach  ntwbar  instead  of  straaa  velocity  corresponding  to  the 
relation: 

P  =  (3-32a) 

'  Ma  d$ 
c 

The  last  tarn  in  itq.  (3-31)  vanishes  again  because  of  the  condition 
A:  V  1.  Fig.  7  has  been  taken  froa  ttie  solutions  obtained  by  Cohan  and 
Rashctko  and  indicates  the  dinensionlass  enthalpy  gradient  at  the  wall 
surface  plotted  over  the  pressure  gradient  paraaeter  p  ,  The  wall 
tei^Mraturv  is  assuaed  constant  in  these  solutions.  It  can  be  c'  served 
to  be  a  characteristic  feature  of  the  diaensionless  enthalpy  gradient 
that  it  varies  only  Moderately  with  the  pressure  gradient  paraneter 
P  .  oapecially  for  small  values  of  the  parameter  the 

variation  become  quite  small.  A  small  value  of  the  parameter 
refers  to  a  condition  when  the  temperature  of  the  surfi.ee  is  small 
as  compared  to  the  temperature  in  the  stream  outside  the  boundary  layer. 
Such  a  condition  is  usually  found  in  high  velocity  engii^eerlng 
applications  since  high  velocity  flow  is  goiwrally  connected  with  high 
temperatures,  a;kl  since,  on  the  other  .hand,  the  vail  temperature 
cannot  exceed  a  limiting  value  determined  by  the  strength  characteristics 
of  the  material.  Consequently,  Lees  proposes  to  neglect  i;be  dependence 


32 


'f  the  dtaensionless  enthulpy  gradient  oa  tne  pressure  ^id.ient  for 
such  coniitions  ani  to  use  an  average  value  equal  0.5. 

The  local  heat  flux  per  unit  area  into  the  wall  surface  can  then 


be  expressed  in  the  following  way: 


JSL _  (  la.) 

I-  - 


<iw  ■  V  ■  r"  ^  ^  I  '  ~- 

^  w  pw  '*  »»  ^  ^e  " 

-  ‘Sr  (r  .  )  ?£  (l.-l^J0.50  x 

"v  f«  V?- 


V3  ?.Vo 


(3-33) 


The  factor  ’***  added  in  the  last  tei»  to  account  for  rrandtl 

nuobers  sowwhat  different  fron  one.  The  enthalpy  difference 

..  will  be  replaced  by  the  difference  between  recovery  and  wall 
enthalpies  Ij.  -  ^  ^0/  the  saae  reason,  for  a  highly  cooled  wall, 
where  the  wall  teaperature  Is  ouch  lower  than  the  reco*.-  .0  ter:;'  .rature, 
this  is  a  sofflelently  good  approxlaation.  With  the  definitions  of  a 


heat  transfer  coefficient, iwiseelt  nuaber,and  a  Reynolds  nuiter: 

.  *  —23-.  8e„  r -  (3-34) 

Wli  I,  Olt  Jl, 

“w  CO 

(i... reference  length,  o  refers  to  conlltlon  at  stagnation  point  oatslds 
the  boundary  layer)  the  following  ecoations  are  obtained  fTon  the  above 


relation  for  the  wall  beat  flux: 


Rtt^sO,3S  f(8)  (rtr)  FHBj^ 


‘e  H 


r(s)  * 


U 


ft  ^  -51^“ 


(3-35) 


(3-36) 
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«i'-h  these  e  juatiins  the  local  heat  flux  along  the  surface  of  a  blunt 
object  can  be  calculated  as  soon  as  the  variation  of  de.^ity,  viscosity, 
and  velocity  outside  the  boundary'  layer  and  along  the  surface  Is  kncvn. 
How  to  obtain  these  values  has  been  discussed  in  a  previous  section. 

Any  ai*bitrary  length  can  be  selected  os  reference  length — for  instance, 
the  largest  dinension  of  the  object. 

The  isethcd  ss  outlined  up  to  uok  has  the  short-coning  that  it.  is 
based  on  Fig.  7  vhich  holds  only  for  a  gas  with  the  prod>ict,  density 
tines  viscosity,  equal  to  a  constant.  Fbr  a  gas  with  other  density  and 
viscosity  .’’t.lations,  errors  nay  be  introduced  as  will  be  denonstrated 
by  calculating  the  heat  transfer  coefficient  for  the  stagnation  point 
region  of  a  blunt  object  with  equation  (3-35}.  The  following  relatione 
hold  for  rotationally  syiOKtric  stagnation  flow: 

a  =  1  ,  Ug  =  ex  ,  =  s 

In  addition,  the  density  and  viscosity  variation  in  the  close  neighbor¬ 
hood  of  a  stagnation  point  can  be  neglected  (  f g  -  ^ 

Equation  (3*36)  CMi  then  easily  be  integrated,  one  the  foHcuisg 


Introducing  this  relation  into  the  above  equation  for  the  Nusselt 
nunher  results  in: 

.  0.™  ^  A 

k  ?o  *  ’*•<’  I  Ao 


The  equation  can  be  re-written: 
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-  0,70(P*-)^^  .  0.70{Pr)^^  \!£»V 

^  fot  Kf,^o  irw 

A  comparison  of  this  last  sxpresslon  with  the  equation  obtained  by  Fay 
and  Riddell  for  stagnatioi  point  Q.ew  and  listed  cn  page  23  indleates 
that  equaUon  (3-35)  will  not  properly  describe  the  heat  tr«»sfer  eoiw. 
dition  for  a  gas  in  which  the  product  dmslty  tiaee  Tlscosity  fcUces 
a  relat.  on  which  is  different  from  a  constant*  Lees  corrected  this  slta» 
atlon  by  multiplying  equation  (V35)  with  Ss.^P,  which  does  not  change 

the  relatitn  for  a  gas  with  co»at.  In  Bef,  k37,  it  is  proposed 
to  aodify  the  nethod  ly  coabining  it  with  the  reference  enthalpy  procechm 
Tne  eolation  (>35)  la  raplaoad  by  the  foUcaring  aquation 


V3 

Wtt  -  0*35  r*{s)  (Pr*) 


(3-37) 


in  ^ch  Nnaselt  and  prandU  aMhen  as  well  aa  F(a)  are  r-sferred  to  a 
referenoa  oondLUon  as  dascribed  by  eqiationa  (3^36)  cr  (>17)  md 


based  on  locsl  conditions  at  tbs  surf aos  point  nndar  eoasidaratlon* 
Tha  function  F*(s)  la  dsflned  as 


It  is  ahowa  in  Rsf*  li37  that  the  proposed  aathod  sgrses  very  wsU  with 
boundary  layer  solutions  «id  eoparlnents  reported  In  Rtf*  186* 


A  fiarther  extension  of  the  astfaod  is  tentatirely  proposed  for  si«>- 
uatieint  t<r  Uhlch  the  ratio  is  not  ssall  as  oowparsd  to  1*  m  this 
a  value  which  replaees  the  constant  0*35  in  aquation  (>37)  can  be  read 
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off  Fig.  7  for  the  appropriate  value  of  the  pre&sure  grad'ent  paraMeter 
P  .  This  later  paroueter  is  obtained  from  equation  (3-32a).  itelations 
are  developed  in  Beferences  213>  125*  and  which  are  even  easier  to 
evaluate  than  equations  (3-36)  and  (3-37).  It  is  to  be  expected* 
however*  that  the  aethcd  presented  in  this  secticn  can  be  extrapolated 
with  nore  confidence  to  situations  apart  froit  the  reentry  problee. 
turbulent  Boundary  Ijiarer. 

Approxinate  netbods  sinilar  to  those  used  for  lamuiar  boundaiy 
layer  f-ow  have  also  been  developed  for  the  condition  where  the  boundary 
layer  has  becone  turbulent.  Ihese  methods  have  been  discussed  and 
recently  coapared  in  References  31t5*  itl9*  and  U35.  Agreement  aimng 
the  various  nethods  is  still  rather  poor.  It  has*  however*  been  found 
that  the  influence  of  the  pressure  gradient  on  heat  transfer  is  less 
in  a  turbulent  boundary  l<grer  than  in  a  laainar  one.  Use  of  equation 
(3-l<-')  describii^  heat  transfer  on  a  flat  plate  gives  a  fairly  food 
approxinati  n  for  a  turbulent  boundary  layer  on  a  blunt  object  when  the 
velocity  u^  appearing  in  the  dtanton  and  Reynolds  numbers  is  interpreted 
as  the  local  value  outside  the  boundary  layer*  and  when  the  distance 
X  in  the  deynolds  number  is  interpreted  as  the  distance  s  of  the 
surface  point  under  the  consideration  from  the  stagnation  point  measured 
in  flow  direction  along  the  body  swface. 

312u.  separated  Flow.  Heat  transfer  in  separated  flow  regions 
like  the  one  sketched  in  Fig.  U  is  still  very  incompletely  vaderstood. 
This  region  a  s  cither  more  or  less  stagnant  or  filled  with  a  regular 
circulating  lamiasr  flow  or  with  a  turbulent  flow  containing  irregular 
vorticity  and  of  an  unsteady*  flnetaating  character,  khlch  of  these 
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flow  situations  occurs  depends  on  tne  .te/nolds  rmmuer.  In  a  constant 
property  fluid,  heat  transfer  to  the  surface  of  an  object  bordering  on 
the  separated  region  may  be  smaller  or  larger  than  the  heat  tran-<^fer 
tc  the  front  portion  which  is  covered  by  a  boundary  layer,  ihe  most 
extensive  information  is  available  for  heat  transfer  in  flow  normal  to 
a  cylisdur  with  circulax-  uross-section.  It  has  been  found  that  for 
.-ieynolos  nnbers  up  to  10  the  local  heat  transfer  is  a  minimum  at  the 
reaveard  stagnation  point  of  the  cylinder.  At  hi  'er  Reynolds  numbers, 
a  maxixaim  appears  at  the  rearward  stagnation  point  which  gradually 
increases,  and  for  Reynolds  numbers  oeyond  approximately  50,000  the 
haat  transfer  coefficient  on  the  rearward  portion  of  the  cylinder  it 
as  large  as  the  one  on  the  front  portion,  with  Ibrther  increase  of 
deynolds  nunber  the  heat  transfer  becomes  larger  in  the  back  portion 
than  in  the  ffMnt  portion,  ror  high  velocity  flow  of  a  compresuihle 
fluid,  heat  transfer  in  the  separated  region  is  influenced  by  another 
factor  which  tends  to  decrease  it  relative  to  the  values  in  an 
incompressible  fluid.  It  is  known  that  heat  transfer  in  forced  con* 
vection  depeiads  mainly  on  the  uroduct  of  density  times  velocit.,  and 
in  high  velocity  flov  the  density  in  separated  regions  is  usually 
consiucrably  smaller  than  the  density  along  the  forward  portion  of 
the  object.  Correspoiidingly,  it  was  found  that  heat  transfer  coefficients 
in  separated  regions  of  objects  exposed  to  sur^ersonic  flow  are  quite 
small  as  compared  tc  heat  transfer  in  the  front  portion.  They  are 
osoally  of  the  sane  order  of  nagnit>ice  os  local  heat  transfer' 
coefficients  which  exist  on  the  object  Just  upstream  of  the  point  of 
separation. 

For  sufficiently  small  Heynolds  numbers  for  which  the  flow  in  the 
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.eoarateci  region  exhibits  a  steady  regular  circtilation  of  laminar 
character,  heat  traasfc.  coefficients  were  recently  obtained  by 
analysis  (aefe'vnces  51,  205).  For  higher  iteynclr^s  numbers  where 

the  flow  exhibits  the  irregularly  fluctuating  turbulent  character, 
heat  transfer  coefficients  are  only  ob-.ainable  up  to  now  by  experinents 
and  a  pear  to  be  very  seniitive  to  the  pcrticular  shape  of  the  object 
which  is  exposed  tc  flow. 

32.  illSaOCIttllci.,  10!»iiU.I.N 

file  discussion  up  to  now  was  concerned  with  a  "sixigle  ccaponent* 
gas;  this  means  with  a  gas  in  which  no  changes  in  chemic  d  con.oosition 
occur.  In  reality,  the  tei^erature  increase  occurring  in  high 
velocity  flow  is  often  so  large  that  the  gas  dissociates.  This  is  especially 
true  in  the  regions  behind  snock  waves  on  objects  like  the  one  shown 
in  Fig.  U.  ibe  surface  tenueratures  of  the  objects  are  generaljy^ 
much  lower,  and  therefore  a  teHoerature  drop  in  the  direction  towards 
the  surface  occurs  in  the  boundary  layer.  Accordingly,  the  dissociated 
atoms  recombine  again  within  the  boundary  layer.  Tc  which  degree  they 
can  do  this  depends  on  the  tiiSi,  which  they  have  vailahle.  It  is 
convenient  to  study  two  limiting  cases  which  must  bracket  the  situations 
as  they  actually  occur,  ^ne  limit  is  a  condition  in  which  the  recoa> 
bination  rates  are  very  fast  compared  with  the  time  which  the  molecules 
or  atoms  need  to  change  location  in  the  boundary  layer  by  diffusion  and 
convection.  Accordingly,  themodynamic  equilibriu*  will  be  established 
at  each  point  and  the  dissociation  within  the  boundary  layer  is 
determined  coo^letely  by  the  local  tes^x-atures  sued  pressures.  Tbe  other 
lindting  situation  is  encountered  when  the  dissociation  rat'.^  ere  ver.-  slow 
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compai^  with  the  difiusion  process.  In  such  a  .-ase  the  aolecules  will 
not  find  sttfficlant  tine  to  recoi-tbine  within  the  boundary  lay^r  itsslf . 
ioat  concentration  field  of  the  atoas  is  established  depends  then  on 
the  cv/ndition  at  the  wall.  nOten  the  wall  is  catalytic^  then  the 
reyi.ubi  natiofi  of  aeoas  occurs  at  tbs  wall  to  such  a  degree  that  sqail- 
lb^llnn  is  again  sstahlishsd  is  thaw  gas  at  that  point,  ihe  recoinbinsd  aole- 
eules  diffuse  away  froa  the  wall.  On  the  other  hand,  the  atoms  haws 
to  diffuse  towards  the  wall,  and  this  diffusion  process  determinss  the 
cunesntration  profile  within  the  boundary  layer.  If,  on  the  other 
iiaod,  the  wall  is  not  catalytic,  then  the  recoid>ination  rate  will  be 
slow  even  at  the  surface,  and  the  atoae  will  aaiutaiu  the  concentration 
which  they  had  outside  tha  boundaiy  layer  through  tne  whole  region  up 
to  the  s>*cfa<^  of  tha  object. 

It  isii^t  be  expected  that  the  beat  transfer  to  the  well  of  an 
object  in  high  speed  flow  is  influenced  strongly  by  thi*  dissociation 
and  r«coid>ination  process.  A  closer  inspection,  however,  will  show 
that  this  is  not  the  case,  for  this  puroose,  ws  will  consider  a  gas 
consisting  of  atoms  and  molecules  of  the  same  species  and  at  first 
determine  the  equation  which  govern  the  beat  flow  in  such  a 
dissociating  and  recoahlning  gas.  In  this  situation,  two  processes 
cause  a  heat  flow:  first,  the  usual  thermal  conduction  process,  and 
secondly  the  diffusion  process  of  the  atoms  and  molecules,  because  in 
this  process  the  particles  cany  along  their  enthalpy,  ibe  equation 
de.'cribing  the  beat  flux  q  in  such  a  fluid  with  a  gradient  of  tenperuture 
and  concentration  in  y  direction  reeds: 


-  .k  ^  -  p  0,  I,  “nr 
*  dy  T  ^  ^ 


(3-36) 
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k  in  this  equation  denotes  the  heat  conductivity  (exactly  the  uarameter 
detemiining  the  trunsi^ort  of  energy  according  to  the  translational^ 
rotational,  end  vibrational  degrees  of  freedoa).  is  the  diffusion 
coefficient  for  interdiffusion  of  the  atoms  and  molecules,  is  the 
energy  of  dissociation  per  tinit  mass  of  atoms,  and  is  the  mass 
fraction  of  the  atoms  in  the  gas.  i‘ho  equation  therefore  is  written 
in  mass  frsetion  rather  than  in  concentration.  The  enthalpy  of  the 
molecule -atom  mixture  is  described  by  the  following  equation: 


di  :  c  dX  +•  i,  iw 

p  a  A 


(3-39^ 


in  wtiich  c^  again  denotes  the  specific  heat  at  constant  pressure 
comprising  the  translational,  rotational,  and  vibrational  degrees  of 
freedom,  however,  exclU'lli^  cbsnlcal  energy.  Introducing  this  relation 
for  the  enthalpy  into  equation  (3>3&i  gives  for  the  heat  flux  tbs 
following  relation: 

This  equation  can  be  simplified  by  introduction  of  a  dimensionless 
pei'ameter  celled  i<ewi8  number  and  defined  by: 


irn 


(3-W>) 


With  the  lewis  oumber  the  equation  for  the  heat  flux  reads: 

The  lewis  number  for  gases  has  a  value  which  is  not  too  far  from 
one.  Therefor",  a  good  approximate  infozmation  on  beat  flow  under  the 
influence  of  diwoclatlon  can  be  obtained  when  the  condition  in  a  gas 
with  a  lewis  number  equal  one  is  considered.  For  such  a  gas,  the  heat 
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flow  aquation  sispllflea  to: 


q 


k  ^ 

■  Cp  'ly 


(3-li2) 


It  c'4a  ba  8««i'2  that  tha  •quatioa  for  tha  liaat  flow  la  Juab  aa  aiapla 
as  for  a  alngla  componant  gaa.  iha  only  dlffaranc-e  is  eausad  by  tha 
fact  that, in  th/j  enthalpy  gradient  daaeriblng  tha  heat  flow,  tha 
chaalcai  energy  (energy  of  diasooiatioiO  baa  to  ba  ineludad.  On  tha 
othar  hand,  it  should  ba  kept  In  mind  that  tha  heat  eonduotiTitqr  and 
specific  heat  do  not  include  tha  ehaaieal  energy  but  only  the  pure 
gas  energies.  One  haa  to  distinguish  those  properties  fron  tha  ones 
inoluding  tha  ehonloal  energy  which  often  are  Hated  in  refereaoe  works* 
▲  more  detailed  discussion  of  this  diffarenea  will  be  glean  in  a 
later  section.  Equation  (3*U2)  indicates  that  tha  heat  flux  la 
dependent  on  the  enthalpy  gradient  in  the  gas  alone*  end  that  It  does 
not  matter  at  all  whether  the  tcanspert  in  detail  is  oaused  bf  oonduc- 
tion  or  by  diffusion.  IMs  already  indicates  that  the  heat  flow  does 
not  depend  oa  the  degree  of  dissociation  in  the  gae.  An  example  will 
explain  more  clearly  what  is  meant  by  this  statement,  let  us  eenelder 
for  this  purpose  a  gae  layer  with  finite  thickness  b  and  let  us 
simplify  the  actual  conditions  by  assuming  that  hsat  conductirity 
and  specific  heats  are  constant.  In  this  case  the  equation  oan 
inmediately  be  integrated  over  the  thiokneea  of  the  layer,  and  the 
following  equation  results: 


_  k 

q  .  ,  ,  - - - 

P  h 

i^  and  i  denote  the  enthalpies  at  the  two  borders  of  tha  layer.  The 
notation  is  used  because  the  layer  nay  be  considered  as  a  crude  model 
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of  a  boundaiy  layer,  then  corresponds  to  the  enthalpy  of  the  gas 
at  the  solid  surface,  and  to  the  enthalpy  at  the  outer  edge  of  tha 
boundary  layei.  Ihis  enthalpy  i^  is  prescribed  through  the  concentra¬ 
tion  of  atoBS  and  through  the  temperature.  The  entlialpy  i^  at  the  wall 
surface,  on  the  other  hand,  depends  on  the  reccnbination  conditions  at 
the  surface,  if  the  wall  is  catalytic,  then  the  concentration  of  atoms 
is  equal  to  the  equilibriuit  concentration  which  belongs  ou  the 
temperature  and  pressure  of  the  wall  surface.  At  a  prescribed  wall 
temperature  and  pressure  in  the  boundaxy  layer,  the  enthalpy  is  also 
a  fixed  value,  and  this  is  so  regardless  of  the  fact  whether  the  atoms 
recombine  within  the  laye»*  itself  or  not  —  in  other  words,  whether 
equilibrium  or  froten  state  or  any  condition  in  betv.een  exists  within 
the  layer.  This  then  indicates  that  tne  heat  flux  to  tr.e  wall  at 
prescribed  wall  temperature  and  prescribed  coixiitions  outside  the 
layer  is  a  fixed  value  and  independent  of  the  ch  ■.  omb'  laticn 

process. 

If,  on  tne  ether  hand,  the  wall  is  non-catalytic,  then  the  heat 
flux  depends  on  the  condition  for  recombinaticn  within  the  layer. 

For  equilibrium  state  the  concentration  within  the  layer  is  everywhere 
equal  to  tne  equilibrium  concentration  belori;;ing  to  the  local  temperature. 
Correspondingly,  the  equilibrium  concentration  is  established  at  the 
wall  surface  by  the  conditions  in  the  gas  Itself,  and  tne  beat  flow 
will  be  again  the  same  for  a  catalytic  wall.  If,  on  the  other  hand, 
the  reccedhination  in  the  gas  is  negligible  (frozen  state),  then  the 
concentrations  will  be  uniform  throughout  the  whole  layer  of  the  gas 
and  will  also  not  be  influenced  by  the  presence  of  the  wall.  In  this 
case,  no  transport  of  enargy  by  dif fusion  occurs,  and  only  the  first 
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tarn  in  equation  (3-38)  is  active  tor  the  transport  of  energ/.  For  a 
constant  heat  conductivity  this  term  can  again  be  intagrateci  ^  gives: 


It  is  clear  that  the  boat  transport  in  this  case  will  be  less  than  when 
diffusion  is  present.  The  ratio  of  the  heat  flux  in  a  frosen  state 
gas  and  at  a  non-catalytic  wall  to  the  heat  flux  at  a  catalytic  wall 
ccn  bo  expressed  in  the  following  way: 


1.  . 


S,  ‘v-i. 

Note  that  c  only  oontain*  the  translational,  rotational  and  vibrational 

P 

energy,  whereas  1  contains  additionally  chemical  energy. 

In  a  boundary  layer,  the  energy  transport  by  conduction  and 
diffusion  is  ;tlll  described  by  equation  (3* 3b).  To  It,  however,  has 
to  be  added  the  transport  by  convection.  Mumaroue  ea*  dilations  have 
been  performed  In  which  the  laminar  boundary  layer  equations  hsve 
been  solved  for  the  determination  of  heat  transfer  in  a  high  velocity 
air  flow  undur  the  preaenoe  of  dlaeoclatlon  end  recosbination*  The 
corresponding  reports  are  listed  In  the  list  of  references.  FUy  and 
Hidden  (Ref,  i2l|)  investigated  the  heat  transfer  In  a  rotatlonally 
syawtrlc  stagnation  flow  of  air  and  arrived  at  the  statement  that  heat 
transfer  at  a  catalytic  surface  Js  obtained  from  the  corresponding 
equation  in  a  single  component  gas  by  aniltlplleatlon  with  the  following 
factor : 


(U“  -  1) 


^•D 
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r,g.s 

Temperature  arsd  concentration  profiles 
on  catalytic  and  noncatalytic  surfaces 
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If  our  relation  v3-3Z^  is  used  then  the  foUoHiag  equation  it.  obtained 
which  describes  heat  transfer  in  rotationally  sywnetric  stagnation  flow 
of  a  dissociating  air  boucdaiy  layer: 


1 


i 


ijal.  :  0.70  (ft-*) 

HP 

Indicates  the  energy  of  dissociation  per  unit  mass  of  the  gas 
consisting  of  molecules  and  stoias,  and  i^  is  the  enthali^  of  the 
sixtura  at  the  stagnation  point  and  outside  the  boundaiy  liy^r.  Ihe 
expoiient  n  was  detemined  to  have  the  value  0.52  for  equilibrium  state 
and  0.63  for  froaen  statf  within  the  boundary  layer.  For  Couotte  flow 
a  value  0.5  had  been  fou»l.  "t  can  therefore  be  assumed  that  the 
exponent  n  depends  only  moderaMiy  on  the  specific  flow  situation  and 
that  the  above  correction  factor  can  generally  be  used  to  obtain  heat 
transfer  in  a  dissociating  or  recombining  gas.  The  Lewis  number  In  air 
for  h/P«:^8onic  flow  conditions  is  expected  to  assume  vfl.  hs  between  1 
and  1.5>  and  a  consideration  of  diffusion  and  recombination  m«y  give 
an  increase  of  the  ineat  transfer  up  to  approximately  20  per  cent.  On 
the  Cbher  band,  for  frosen  state  in  the  boundary  layer  azxi  a  non- 
catalytic  wall,  a  possible  redaction  of  the  beat  transfer  parameter 
to  one-third  of  the  value  which  it  would  have  if  the  wall  were 
catalytic  was  calcvQated  by  Fv  and  Hiddell  iHef.  'Xhia 

calculation  was  mode  for  a  iligbt  condition  as  It  nay  occur  in  the 
re-entry  of  a  missile.  Frosen  state  in  the  boundary  layer  has  been 
essumed,  and  this  requires  that  the  flight  occurs  at  very  high  altitude.' 
Fig.  8  gives  a  schematic  sketch  of  the  variation  of  temperature  and 
concentration  through  a  boundary  layer;  once  for  a  catalytic  and 
second  for  a  non-catalytlc  surface.  The  temiwrature  profile  In  the 
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bouiKlaiy  layer  is  essentially  unchanged  by  the  surface  condition.  On 
the  other  hand,  the  concentration  is  uniform  for  a  non-catalytic 
surface  and  decreases  to  its  equilibrium  value  according  to  the  lower 
surface  tempereture  for  the  catalytic  surface,  energy  transport  is 
affected  by  the  temperature  gradient  alone  at  tte  non-catalytic  surface 
and  by  tiie  concentration  gradient  and  temperature  gradient  at  the 
catalytic  surface. 


U.  1U:>S  ftUajyiffl  COOLIHQ 

Ibe  general  equations  which  describe  the  tamnerature  condition  in 
a  solid  wall  and  the  heat  Ilux  to  the  jurface  of  such  a  wall  under  the 
influence  of  mass  transfer  cooling  have  been  discussed  in  cliapter  2, 

In  the  present  chapter  we  are  coneamed  only  with  the  process  within 
the  boundaiy  layer  itself,  which  determines  the  heat  transfer  coefficient. 
Correspondingly,  the  term  mass  transfer  cooling  will  be  interpi'eved 
here  as  ir.compi9slng  all  cooling  methods  by  which  a  mass  flow  away  from 
the  surface  and  into  the  gas  stream  is  generated,  regardless  of  the 
specific  method  by  which  such  a  mass  flow  is  achieved.  It  includes 
processes  like  transpiration  cooling,  sweat  cooling,  ablation  cooling, 
liquid  film  cooling,  and  so  on.*^  The  coolant  which  in  gaseous  form 
moves  twsy  from' the  surface  of  tl»e  solid  object  is  usually  of  a 
substance  different  from  tiia  fluid  aovii^  in  the  outside  flow  over  the 
surface.  Correspondingly,  the  mass  transfer  cooling  process  will,  in 

The  interphase  between  liquid  and  gas  in  film  cooling  wves 
along  the  suv-faoe.  The  welocities  of  this  shear  flow  in  tbs 
film  are,  hovever,  so  small  that  their  effect  on  the  gas  boundary 
layer  Is  negHgible. 
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addition  to  ll«t  paraaeters  occurring  on  a  aclid  '’urfac*,  depend  on  tvo 
new  psraneters  *■  on  the  enount  ot  ness  released  at  the  surrace  end  on 
the  nature  cf  .ne  cool&nt  gas.  eccordlngljr,  it  requires  a  nueh  larger 
nunber  of  ralculatluns  or  experlaents  to  obtain  the  required  knoi^edge 
on  the  Diasa  transfer  cooling  process.  This  Is  the  reason  idQr  our 
present  day  infontation  on  nass  transfer  coollnc  is  still  aifvne  ss 
coi^>ared  with  that  on  beat  transfer  on  a  solid  surface.  Osoeralljr,  it 
is  found  that  a  nass  flow  away  fron  tba  surface  raduees  heat  transfer. 

Ibe  procedure  which  appears  bast  in  the  light  of  tbs  renartes  nada  aboes 
and  which  today  has  been  quite  generelly  eecepted  Is  to  calculate  hast 
transfer  on  a  uass  cooled  wall  In  tbs  following  way:  Ibe  heat  transfer 
on  a  solid  surface  under  equivalent  flow  conditions  is  calculated  and 
this  heat  tran£far  coeffieiant  is  then  nultiplisd  with  a  eorrsetlOB 
factor  tduLCh  describes  the  reduction  obtained  tqr  the  ness  transfer 
cooling  process.  Ibis  procedure  will  be  need  in  the  foUoMing  dieeueslon. 

la.  WITHOUT  CH&ilCAL  RKACXIOIIS 

In  Many  apolications»  the  cheaical  speeiee  released  at  the  surfeee 
will  react  with  the  conponents  in  the  outside  gss  flow.  Ibe  dlsensslcm 
in  this  section  will  assuns  that  such  resctloas  do  not  occur.  In 
addition,  it  will  consider,'  on  the  one  hand,  the  gas  in  the  outside 
How  and,  on  the  other  band,  tha  nass  released  tnm  the  surface  aa 
single  conponents*  This  is  adnisslble  as  long  a«  each  of  tha  tm 
conponents  consists  of  species  whose  properties  differ  ew^Mratively 
little.  'Ws  are  therefore  concerned  with  the  nass  aovansat  bgr  oonrsetioa 
and  diffusicn  of  two  coq^onants  in  a  gas  nixtux*  r^tiv*  to  each  other. 
It  will  be  seen  that  this  process  has  often  a  prcaocccad  affect  on 
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iieat  transfer.  The  discussion  viU.  again  be  earned  out  in  the 
following  order,  at  first  a  Mixture  of  two  cuapoosnts  will  be  conhioered 
with  pNOperties  which  are  constants  (indepenoent  of  pressure, 
teaperatur«>  and  coaposition.).  as  a  next  step,  the  prope:*tie8  of  the 
Mixture  will  be  considered  to  depend  on  tenpsrature  but  not-  on 
CuBpositlun.  Such  a  situation  is  q-U.te  often  approxhaated  in  actual 
Mass  transfer  cooling  processes.  In  a  third  step,  the  properties  will 
be  considered  to  depend  on  t^paratur*  srd  ca  coaposition,  ami  in 
addition  -o  be  widely  different  between  the  tve  cosponants.  In 
recent  investigations,  it  has  be«)n  fouitd  that  gases  with  low  Molecular 
weight  appear  to  offer  special  advantages  as  coolants  for  a  itass 
transfer  cqpling  process,  jince  the  properties  denend  essentially  on 
the  Molecular  weight,  such  a  situation  has  created  interest  in  the 
influence  of  large  differences  in  properties  between  the  coolant  and 
the  outside  gas  on  the  mass  transfer  coolii^  process. 

lill.  Constant  ftoperties 

The  boundary  layer  equations  for  flow  over  a  surface  with  mass 
release. _are  .the  sane  as  those  describing  the  velocity  and  tenperature 
fields  on  e  solid  wall.  The  only  change  occurs  in  the  boundaiy 
conditions  where  a  finite  velocity  v  -  m/^  nonaal  to  the  surface  is 

now  prescribed.  Solutions  of  these  equations  for  laeinar  and  turbulent 

^  • 

boundary  layei'S  with  oass  transfer  cooling  hare  shown  that  Stanton 
nuaber  or  liusselt  nuidwr  describing  the  local  heat  transfer  coefficient 
depend  on  the  following  paraMeter 

(Hs)“  ik-li 

V"e 
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when  the  mass  m  released  per  unit  time  and  surface  area  raries  ioversel/ 

It 

prcportional  to  the  expression  x  u^'  .  Ihe  exponent  n  has  the  value 

U.5  for  a  laminar  boundary  laymr  and  0.?  for  a  turbulent  one.  For  this 
situation,  the  boundary  layer  profiles  at  different  locations  x  ur 
similar,  oset  transfer  coefficients  obtained  from  these  sladlerlty 
solutions  are  often  used  as  aoproximatlons  for  a  determination  of  heat 
transfer  amen  wben  the  mass  flow  distribution  omor  tho  surfsce  is 
different  fm*  the  css  cn  which  tbs  calc^aiion  was  based.  Ibis 
procedure  Is  snelototta  to  tbs  nstbod  of  "loesl  slatlerlty*  whl«^  hss 
boon  sppliod  In  Saetlon  5123.  For  taibulsat  flow,  it  can  be  ei^eetsd 
tbat.  the  influsoeo  of  e  speelfic  msss  relesse  distribution  on  heat 
transfer  is  eomparatlmely  smell.  For  laminar  bouadaxy  Isyers  the 
iaflueace  may  be  consldarable.  Itomeer,  little  lafoxaetlon  is  evell* 
aUle  from  which  its  ■egiitnrte  can  bo  estimated. 

fbe  parameters  may  be  brou^t  to  a  fozm  which  Is  ideatlcel  for 
i—<«er  aarl  turbulent  flow  wten  the  dqraolde  namber  is  replaced  by 
Stantons  nueber 


(U-2) 


lbs  Stanton  rustier  St^  deeerlblac  beat  tresefer  on  a  solid  surfaoe 
under  the  asms  bonadsry  eondltloas  (eme'<pt  ▼  S  0)  le  used  fbr  this 
peimmeter  in  the  present  report.  C.  •).  Seott  demonstrated  additlonelly 


that  the  results  of  eaZouletioae  on  beet  transfer  in  leminer  boundary 


layers  with  preseare  (radlente  can  be  brou^it  to  near  eoincidende 
wbra  the  psrnMter  (U-2)  is  ehsacad  to 


WADC  TR  59-62U 


49 


I 


V 


Reducti''-  of  hnat  trnn’^l’er  and  rooovery  fac*cr  by 
tranefer  co-'llfiff  (or.stnnt  oropcrV  fluid) 


is  th*  pressure  gradient  paramter  a«  useu  in  SecUon  3123.  It  bw 
the  value  0  for  flat  plate  flow,  1  for  plane  stagnation  flow,  aal  0.5 
for  rotationilljr  sjnwetric  ategnaUon  flow,  i'or  turbulent  flow, 

/3  -  0  sQould  be  a  good  approximation  even  when  pressure  gradients  are 
present. 

Figure  0  presents  ths  rstis  it/h^  of  the  actual  heat  transfer 
coeffieiant  to  the  one  on  a  solid  surface  under  the  sim  outside  flow 
conditions  as  a  fraction  of  the  nesa  relaaee  ner«ee»jM»  (|;.3)  ^ 

^  ~  0,?,  It  aaF  be  observed  that  the  values  for  flat  i:d.4te 
flow  and  for  plane  stagnation  flow  agree  quite  eloseljr  in  this  present¬ 
ation,  The  figure  also  contains  ths  ratio  r/r^  of  the  eetual  teeperetnre 
recovery  factor  to  the  one  on  a  solid  surfhc#  for  1—inar  flat  plate 
flow.  Knowledge  of  the  recovery  factor  is  of  nlnor  Inporiance  In 
■ass  transfer  cooling  because  the  difference  betwaen  total  gas  twi isrituis 
and  i«ll  te^rature  Is  usually  large  ecapared  with  t'e  difference 
between  total  gas  teaperature  and  recovery  teaperature,  an  error  in 
the  recovery  factor  then  has  litUe  influence  on  the  calculation  of 
the  heat  flux. 


iil2.  i^pertief  Deperrdent  on  i^»^Tature 
It  has  been  deaocstreted  in  Reference  lr25  that  the  results  of 
nmerous  leainar  boundary  Uyer  solutions  for  flow  over  a  flat  plate 
•free  with  the  constant  property  roletlon  presented  in  ilg.  9,  wt'tn 


the  blowing  paraater  is  changed  to  the  following  fom  euggasted  by 
H.  Saron: 
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This  paraiaeter  can  be  tratisionned  in  the  following  way: 


This  indlcatas  that  toa  rule  proven  so  successful  for  heat  transfer  on 
solid  surfaces,  nanejor,  ihai  properties  should  be  introduced  ot  a 
reference  taaperature  applie.  to  the  nass  release  paraaeter  as  wall, 
the  oaraaeter  can  again  be  generalited  by  intrcducticn  of  the  Stanton 
nuhber  and  by  addition  of  the  pressure  gradiant  paraaeter.  It  can  be 
espectcd  that  the  paraaet?r 


correlates  heat  transfer  coefficients  for  laalaar  boundary  layer* 
with  pressure  gradients  and  for  turbulsat  booadaiy  layers  (setiinsj)-  0). 

9  can  than  be  used  to  calculate  heat  transfer  in  fluid  with 
taaperature  dependent  nr  parties  idisn  ths  properties  in  ths  aass 
relsast  paraaeter  are  introduced  at  refereoce  teaiaratare  or  reference 
enthalpy.  The  reference  state  is  calculated  froa  equation  (3>13)  or 
(3-16). 


Fig,  10,  taken  Area  leferenee  U26,  presents  ths  ratio  of  ths 
•otual  hsat  flux  q  at  a  aass  transfer  cooled  wall  to  the  heat  flua 
q^  which  the  solid  well  would  have  under  the  saae  flow  conditions,  ^tted 
over  the  aass  release  paraMter*  It  has  bsen  found  that  ths  results 
of  various  calculations  oonelate  on  those  enrves  with  reasonable 
accuracy.  The  notation  on  the  varir'is  curves  iadieatas  the  nature  of 
the  two  ccaponent  aixture.  The  first  tens  gives  ordati-  gas,  and 
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it  cut  be  obserred  that  the  outside  flow  orer  the  surface  was  air  in 
all  oases.  The  Tarlous  curres  represent  with  good  accuracy  the  results 
of  boundary  layer  calculations  on  a  flat  plate  with  a  large  range  of 
Mach  nua^rs,  wall  temperatures  and  free  stream  temperatures.  It  can 
be  obaerved  that  the  curres  arrange  thamselres  according  to  the 
molecular  weight  of  the  coolant,  in  this  way  it  is  suggested  to 
consider  the  molecular  weight  of  the  coolant  as  the  main  parameter 
detezsiniag  the  coolant  effect  in  the  maee  transfer  process.  In 
Refersnc<-  1x26  *  it  has  been  found  that  the  rarious  eurree  in  Fig.  10 
can  be  expressed  by  the  following  equation: 

q/«j,  a  1-2.1,6(  ?^V?eV  >^a 


* 

^  and  ^  denote  the  deneity  and  riseoBity  of  air  at  rwrerease 
teK)eratare*  in  this  equation  indioatea  the  molecular  weight  of 
air  and  the  molecular  weight  of  the  coolant  gas.  Tht  scm  equswion 
deeerlbee  the  ratio  of  the  actual  wall  ahear  on  the  trenepiration 
cooled  wall  to  the  shear  which  would  exist  on  a  solid  surface  under  the 
esM  flow  condition  iMmb  tte  constant  is  changed  to  a  value  of  2.06. 
a  similar  relation  was  obtained  in  deference  366  for  the  ratio  of  beat 
transfer  coefi'lelents  in  a  laminar  and  torbulsnt  boundary  layer  on  a  flat 
plate.  Tbs  exponent  on  the  molecular  weight  ratio  in  this  reference 
was  found  to  be  0.U  for  laminar  flow  and  0.6  for  turbulent  flow.  Reference 
425  gave  2/3  for  the  exponent  1/3  In  turbulent  flow.  This  result  and  the 
generalised  parameter  can  now  be  used  to  obtain  the  relation 


=  1  -  C  V~  )  (U-7) 

?  ®e  St,  "* 

C  has  the  value  0.73  for  laminar  and  0.37  for  turbulent  flow. 
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The  exponent  n  is  1/3  for  laaiinar  and  2/3  for  turbulent  flow, 
rrooertie-i  for  air  at  reference  temperature  hare  to  be  Introduced  for 
^  *  a«i  in  ^t^*.  The  ratio  Iv'hj  of  beat  transfer  coefficients  is  not 
quite  identical  to  the  ratio  q/q^  of  heat  fluxes  because  the  recorezT' 
temperature  on  a  solid  wall  is  somewhat  diffei-ent  from  The  recofery 
temperature  on  a  mass  transfer  cooled  wall.  In  the  boundary  layer 
solutions  on  which  equation  is  based,  this  difference  is  very 

snail.  This  has  already  been  pointed  out  at  the  end  of  the  pteceedlng 
section.  Reference  1^26  contains  iirfomation  in  fig.  U  on  the 
recovery  factors  of  transpiration  cooled  surfaces. 

Equation  (U-?)  i>  a  relation  which  represents  results  of 
atralyses  and  experinents  available  today  with  an  accuracy  which  is 
sufficient  for  engineering  pur^ses.  However,  it  should  be  rsMidiersd 
that  the  range  of  infomation  to  date  is  still  restricted,  and  that 
con'ections  to  this  relation  aay  have  to  be  made  when  more  information 
becomes  available  in  the  future.  There  are,  for  instance,  som 
Indications  from  calculations  presently  peirformed  at  the  Naval  Ordrvutee 
Laboratory  that  the  molecular  weight  is  actually  ix)t  the  only  perameter 
to  describe  the  behavior  of  a  specific  coolant.  It  may  well  be  that 
in  the  future  either  other  properties  have  to  be  considered  in  the 
relation  (U-7^  or  that  a  procedure  will  have  to  be  used  which  considers 
all  the  properties  of  the  coolant  gas  as  well  up  of  the  gas  in  the  out> 
side  flow  as,  for  instance,  a  reference  concentration  procedure.  Such 
a  reference  concentration  procedure  hoc  bsaa  pro7o«od  in  Reference  li36, 
it  should  also  be  pointed  out  that  equation  (U-7)  is  based  on  analyses 
and  experiaents  in  all  of  which  the  surface  was  inpermeable  to  the 
one  coi^}onent  of  the  gas  mixture,  namely,  to  the  in  the  outside 
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flow.  A  different  sicuation  sometimes  encountered,  for  instaiKe, 

V»^**  A  MW 

the  ojQTSen  neceesary  for  the  reaction  has  tc  be  transported  through 
the  boundary  layer  toward  the  interface  on  which  the  reaction  occurs, 
and  the  relation  between  the  »ass  flow  pcraaeter  and  the  local 
cozKentrations  of  the  two  coapcnents  will  be  different  frcm  the  one 
considered  in  the  prewious  analyses.  This  aay  also  have  an  effect  on 
the  reduction  of  the  heat  transfer  coefficient. 

tt2.  WITH  CUbaCai  REACTIONS 

A  Mass  transfer  cooling  process  is  often  connected  with  chemical 
reactions  when  it  occurs  at  high  temperatures.  These  reactions  change 
the  temperature  field  and  the  properties,  and  Influence  the  heat 
transfer  coefficient  and  the  heat  flux  to  the  surface.  Thi.  -nix^snee 
of  cheaieal  reactions  on  heat  transfer  will  at  first  be  discussed  on  a 
simplified  model  which  was  considered  in  Reference  115  and  l6o.  In 
this  reference,  heat  flow  to  a  surface  in  flat  plate  and  in  plane 
stagnation  flow  was  considered  when  a  steady,  two-dimensional,  laminar 
boundaxy  layer  exists,  and  when  the  flow  Telocitles  arc  sufficiently 
low  to  make  dissipation  negligible.  It  was  assumed  that  the  surface 
releases  a  mass  m  and  that  this  mass  (for  instance,  carbon  or  hydrogen) 
reacts  with  the  ojqrgen  of  the  air  stream  noring  orer  the  surface.  The 
fluid  properties  are  all  postulated  constant  and  hawing  practically  the 
same  walue  for  all  ccmponents  of  which  the  air  and  the  combustible 
material  exists.  Schmidt  number  and  Arendtl  number  hanre  the  value  0.? 
idiieht  means  a  lewis  randiar  equal  to  one.  The  chemicsl  resctisn  'ntea 
are  assumed  to  occur  very  fast  so  thr.t  the  reaction  prccess  is 
comfletely  diffusion  controlled.  Cheidcsl  equilibrium  may  be  close 
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Enthalpy  profiles  in  a  bouivlary  layer 
with  ooafaustion 


■'AhC  TB  ??-A2l4 


57 


to  complete  combustion.  Fig.  11  shous  as  an  example  of  the  results 
of  this  calculation  the  enthalpy  profiles  plotted  over  the  wal?. 
distance  y.  Ihe  wall  temperature  is  assumed  prescribed.  Ihe  full 
lines  irxiicate  the  enthalpy  profiles  for  various  mass  release  rates. 
The  peak  in  the  enthalpy  curves  coincides  with  the  location  of  the 
flame  front  within  which  the  conbustion  occurs.  For  very  small  mass 
release  rates,  the  combustion  occurs  right  at  the  surface,  and  the 
enthalpy  profile  is  the  lowest  one  indicated  in  the  figure.  Beyond 
a  fairly  small  release  rata,  the  combustion  front  lifts  off  the 
surface  and  moves  with  increasing  release  rate  more  and  more  into  the 
interior  of  the  boundary  layer.  The  enthalpies  presented  by  the  full 
lines  constitute  the  sum  of  the  sensible  heats  of  all  the  comoonents 
of  idilch  t^ie  gas  is  composed,  if  the  chemical  nnergy  which  can  be 
released  by  combustion  of  the  local  (uqrgen  is  added,  then  the  dashed 
curves  are  obtained  idtlch  indicate  the  total  enthalpy  profiler  witoLu 
the  boundary  layer.  These  total  enthalpy  profiles  hare  exactly  the 
same  shape  as  the  enthalpy  profiles  for  •.  constant  property  fluid 
without  combustion  and  with  the  proper  mass  release  rate  provided  the 
enthalpy  at  the  wall  is  equal  i^  and  the  enthalpy  at  the  outer  edge 
of  tne  boundary  layer  is  equal  to  the  total  enthalpy  at  that  location. 
Sensible  heat  and  chemical  energy  of  possible  reactions  have  to  be 
Included  in  the  latter  value.  Fr<xn  this  statement  it  becomes  obvious 
that  the  combustion  process  occurring  in  laminar  boundary  layer  flow 
of  a  fluid  with  constant  oroperties  and  L*  s  1  is  included  in  Sts 
effect  on  heat  transfer  when  the  heat  flux  is  calculated  with 
equation  (2-3)  in  which  the  enthalpy  at  the  outer  edge  of  the  boundary 
layer  eomsrises  chemical  energy  as  veil  as  sensible  heat.  The  he;.'. 
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tranafor  co-ifficlent  can  bw  UaUruined  from  «qu»'  ion»  (3-6J  and  (U-7), 
This  procedure  apt'lies  to  a  situation  in  which  the  oxygon  content  at 
the  surface  is  aero.  When  the  oxygen  content  at  the  surface  is 
difiererit  from  zeroi  then  the  chemical  energy  corresponding  to  this 
oxygen  content  must  be  included  in  the  enthalpy  of  the  gas  at  the 
wall  rarface. 

Ihe  explanation  for  this  behavior  can  be  obtained  by  a  reconsider¬ 
ation  of  Section  32  and  specifically  of  the  doTelopoient  which  led  to 
equation  C3-U2).  The  discussion  in  that  section  was  concerned  with 
dissociation.  It  may,  however^  be  applied  in  the  same  way  to  any 
chemical  reaction  within  the  bovundary  layer,  liquation  (3-U2)  than 
Indicates  that  the  heat  flux  normal  to  the  stream  lines  at  ary  point 
within  the  boundary  layer  Is  determined  by  the  gradient  of  the  total 
enthalpy  comprising  sensible  heat,  kinetic  energy,  and  chemical 
energy,  and  that  it  is  imnaterlal  for  this  transport  whalhar  it  ooours 
as  conduction  of  sensible  heat  or  as  diffusion  of  ehemioal  enthalpy. 

In  Fig.  11,  the  transport  of  total  enthalpy  is  purely  by  conduction 
along  the  enthalpy  eunres  in  the  region  between  the  flame  front  and 
the  wall;  it  ia  by  diffusion  as  well  as  by  eonduotlon  along  the  dashed 
part  of  the  enthalpy  curres  in  the  outer  region  of  the  boundary  layer 
outside  the  flans  front* 

L.  Lees  developed  in  Reference  U32  the  laminar  boundary  layer 
equations  for  plans  and  rotationally  synnetrle  flow  of  a  two  component 
gi:i8  mixtuio  with  ohenicel  reactions  and  with  mass  transfer  at  the  well* 
The  continuity  end  momentum  equations  ai^  the  same  as  equations  (3-25) 
and  (3-2<6)*  The  energy  equation  is: 
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•^[r‘^- fe>*‘i4rJ  (W) 

Tbf:  total  anthalpy  1  now  contains  ehosical  energy  In  addition  to  sensible 
beat  and  kinetic  energy.  0  is  the  diffusion  coefficient,  tbs  nass 
fraction  of  the  1^  coaponent  (i  running  from  1  to  2),  is  the 
enthalpy  of  the  i^  conponent.  The  last  tern  in  the  equation  dlsanoears 
when  the  lewis  nuaber  is  equal  one.  The  fact  that  equations  (3-2$)* 
(3>a6),  and  (U-8)  are  not  different  there  fnon  those  for  a  gas  Mixture 
without  cheaical  reactions,  indicates  that  the  velocity  field  and  the 
total  enthalpy  field  are  not  influenced  by  any  oeensring  reactions. 

This  holds  for  a  Mlxtaro  with  constant  properties  and  also  for  a  gas 
with  “  constant  and  constant  ^  when  the  pressure  along  tbo  surface 
is  constant.  The  last  statenent  is  easily  verified  by  inspection  of 
the  equations  (3-30)  and  (3-31)  which  are  the  transfoxaed  aquations 
(3-2$),  (3-26),  and  (3-27).  t')r  s  constant,  C  has  the  value  1, 

**id  for  constant  pressure  ^  is  equal  to  0-  The  aquations  are  then 
the  equations  for  a  constant  property  fluid.  The  enthalpy  gradient  at 
the  wall  surface  deterainss  the  heat  flux  into  the  wall  accordihg  to 
equation  (3-U2). 

The  following  rule  for  a  calculation  of  beat  uransfer  in  aass 
transfer  cooling  and  with  chenical  reactions  follows  froM  these 
considerations:  One  calculates  the  beat  transfer  coefficient  frofai 
the  relations  in  the  preceding  section,  for  instance,  with  equation 
(ii-7).  The  heat  flux  to  the  surface  is  then  calculated  using  the 
right  hand  expression  in  equation  (2-3)  and  Interprstini  enthalpies 
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as  total  eatiuJoles  eoataiplng  seusible  hast,  kinetic  eaarfsr*  mA 
ctwnileal  energy  corresponding  to  possible  xvsctions  within  tbs 
boundsiy  Isgrer.  This  proesdur*  bolds  strictly  for  gases  with  eonstsat 
properties  or  with  ~  constant,  Pr  =  1,  end  constant  pressure  and 
in  anjr  ease  for  le  t  1.  It  should  be  s  good  epproxiastion  when  Br 
sn*  le  ere  not  too  different  frcst  one  and  for  Moderate  pressitre 
gradients,  especially  tdion  proparties  are  again  introdae44  at 
refarenoa  entnalpy. 

The  Msss  releasa  rata  at  tdtieh  ia  required  for  this  eeleolation 
has  to  be  obtained  froA  OTarall  balancas.  In  an  eviration  or 
subliMtion  process,  for  inatance,  it  is  find  bgr  the  energy  helsnee 
in  Chapter  2  together  with  the  tfaeiaodTasHic  relation  betwaaa  avaporeticn 
or  subUaatlon  taaperatnre  and  the  pressure  or  partial  prcsaors  of  the 
released  substance  at  the  surface. 


low  DQidITl£S 


In  the  discussion  up  to  now,  the  fluid  iavclTed  jiut  considered  to 
be  a  eontinouM.  In  reality,  gases  consist  of  indieidosl  adeculss 
and  this  structure  Makes  itself  felt  at  low  densities  which,  for 
instance,  are  obtained  in  the  fli^  of  en  aircraft  at  eery  bi^ 
altitude.  The  paraneter  which  ceteraines  whether  the  Molecular 


structure  influences  heat  transfer  and  friction  is  the  Xnudsen  !iad>er> 


(5-1) 


X  indicates  the  Mean  Molecular  path  length  and  L  denotes  a  eharactar- 


istic  diaension  of  the  object  involved ,  The  flow  in  a  boundary  layer. 
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for  Instance,  can  be  considered  with  good  accuracy  to  be  continuiai 
flow  when  the  ratio  of  mean  molecular  path  length  to  boundary  layer 
thickness  has  a  value  idiich  is  small  as  compared  to  one.  Vidth 
increasing  Knudsen  nuid>er,  the  molecular  structure  makes  itself  felt 
at  first  in  the  inmedlate  neighborhood  of  the  surface  of  the  object 
by  the  fact  that  the  flow  slips  over  the  surface  and  that  the 
temperature  profile  in  the  boundary  layer  has  a  very  sudden  change  near 
the  wall.  This  regime  is  referred  to  as  slip  flow  regime.  Sinultan- 
TOUSly,  an  Interaction  between  shock  and  boundary  layer  usually  begins 
to  influence  the  flow  and  heat  transfer,  tn  the  other  end  of  the 
scele  of  Knudsen  c'wdwrs  (lerge  values;  one  encounters  conditione 
where  the  mean  molecular  path  laogth  is  large  compared  to  any  body 
dimension,  ihis  swans  that,  after  buing  reflected  from  the  surface, 
an  approaching  swlecule  has  practically  no  chance  to  be  reflected 
back  and  to  hit  the  surface  a  second  tine.  Ibis  regime  is  .'cf erred 
to  as  free  molecular  flew  regime.  In  between  this  regime  and  the 
slip  flew  regime  Is  the  transition  regime  in  which  the  mean  molecular 
path  length  ana  the  characteristic  body  dimension  are  of  the  same 
order  of  swgnitude.  The  Khudsen  numb'  r  is  uniquely  related  to  the 
Reynolds  and  Mach  numbers.  As  a  c  •  jquenee,  the  various  flow  regimss 
can  be  indicated  in  a  diagram  in  which  the  Reynolds  nzidwr  is  used 
as  abscissa  and  the  Mach  number  as  ordinate.  Such  e  diagi'am  is 
presented  in  fig.  12  taken  from  Reference  h22.  The  transition  lines 
betvsen  the  regimes  are  to  a  certain  degree  arbltraxy.  Compared  with 
a  similar  diagram  in  Rafermnce  110,  the  limits  between  contlnoun  flow 
end  transition  flow  have  been  shifted  to  the  left,  because  recent 
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exporiMsnta  in  lou  danalty  wind  tunnels  indicated  that  the  contiiviM 
regijw  extends  to  smaller  iiiegraolds  nuahers  than  orlginalij  expected. 

It  has  also  been  pointed  out  in  Aeference  l6h  that  the  proper 
conditions  existing  within  the  boundary  layer  hare  to  be  introduced 
into  iteynolds  and  iteeh  nuaber  in  a  detejaination  whether  the  flow  is 
in  the  ecntinuun  or  slip  flow  roflM.  The  best  Talues  to  use  in  such 
a  case  are  probably  velocity  outside  the  boundary  layer  and  the  properties 
at  the  presrure,  tenperataret  and  concentration  as  they  exist  within 
the  boundary  layer  at  the  surface  of  the  object.  It  has  been  found 
in  Deference  168  that,  in  this  way,  practically  all  conditions  under 
which  cooling  of  nissiles  ai:d  satellitos  beeoaes  critical  are  located 
in  the  continuun  regim.  This  is  very  advantageous  because  heat  transfer 
relations  are  nuch  better  established  for  this  regiaa  than  for  the 
sl^  or  transition  regiae. 

«  very  siagile  relation  for  the  heat  transfer  is  obtained  in  the 
free  noleeular  flow  regiae  for  an  object  flying  with  a  hypersonic 
velocity.  For  such  a  condition,  tbs  aaan  aolacular  velocity  is  saall 
ooapared  to  the  vehicle  speed  V,  and  the  aolecules  can  be  coiuidered 
a^  practically  at  rest.  This  allows  calculation  in  a  aii^a  way  of 
tha  snaber  of  aolaeulaa  which  strike  Uw  vehicle  as  it  aovea  through 
the  gas.  If  Izdioates  the  undisfairbed  density  of  the  gas,  and 
M  its  aolecolar  weight,  then  the  nuaber  of  aolecules  per  unit  voluae 
is  ^m/H.  'iba  nuaber  of  aolecules  hitting  the  object  per  unit  tiae 
is  than  given  by  the  following  relation: 

in  which  a  indicates  the  area  obtained  when  the  vehicle  is  projected 
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onto  c  plane  nonaal  to  xim  lllght  direction,  dech  Bolecul*  etri^cH^ 

the  surface  t.ansfers  the  following  energy  to  tha  object: 

2 


The  tora  a  in  this  expression  tiuicates  t#hat  fraction  of  the 
kinetic  energy  the  aoleeule  transfers.  Ibis  term  is  called  accoaodation 
coef^ciegt.  Ihe  maerlcal  values  for  the  accomodation  coefficient 
have  to  be  determined  by  measurements,  ibey  depend  on  the  nature  of 
thu  gas,  on  the  nature  and  condition  of  the  surface,  end  probably  eleo 
on  the  velocity  V.  Only  limited  information  is  available.  Skperiaen'* 
carried  out  up  to  ■vjw  resulted  in  values  betmseo  0.3  and  1  for  air 
and  for  various'  surface  meterlels.  From  the  two  equations  sUted  abovm 
the  beet  flow  per  unit  erea  of  projected  surface  end  per  unit  time  is: 

ee 

‘*wc  “  f"*  T  (5-2) 

It  is  quite  interesting  to  transform  this  equation  into  e 

number.  Ibe  total  enthalpy  of  the  gee  is  precticelly  only 

energy.  Iberefore,  I  2  ¥^2.  Ihe  enthalpy  of  the  surface  can  be 

considered  as  ve^  aMll  compered  to  this  value,  so  that  the  enthalpy 

difference  is  also  V^/2  and  the  heat  transfer  coefficient  based  on 
enthalplee  Is: 

hjS 

The  Stanton  nmiber  is  therefore  simply: 

=•  (5-3) 

Urn  Stanton  nusber  is  therefore  simply  equal  to  the  acconodation 
coefficient  if  it  is  based  on  conditions  in  the  undisturbed  fluid  and 
on  projectod  area.  Stanton  combers  for  continuum  flow  are  usually 
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bjr  two  or  thriM  orders  of  sagnitode  unaller  tbsa  this  Talus,  this 
iodicstss  the  fact  that  a  boundarr  lagrar  sorrss  as  a  kind  of  insulation 
and  decreases  the  energjr  transfer  between  the  aoleeules  in  the  free 
strew  and  the  object  idten  it  builds  up  around  an  object. 

6.  TBiMSITIOh  TO  TURBDLEMCE 

A  prediction  of  the  location  on  a  bodr  where  the  flow  within  the 
boondarr  layer  changes  fron  lanlnar  to  tarbuleot  state  still  has  to 
be  based  orinarily  on  enoirieal  infomation  obtained  in  wind  tnanel 
experiMOts  and  in  free  flignt  tests,  in  extensiwe  literature  exists 
according  to  the  inportance  of  the  knokledge  of  the  transition  point* 

Ibe  jreader  has  to  be  referred  to  the  literataret  and  only  a  few  i leiiilri 
can  be  nade  in  this  report. 

Stability  theory  predicted  that  the  transition  point  would  be 
dolayed  by  a  drereaae  of  tha  aurfaee  tenperature  so  that  Tsry  large 
transition  ttoynolda  nwbers  are  expected  on  objects  with  cooled  eurfeces. 
A  figure  Indicating  the  results  of  such  calculations  is  containad  in 
Beferenre  110.  These  predictions  hare  generally  been  confliaed  on 

C 

slender  objects.  Transition  dsynolds  txaibsrs  up  to  10  hsre  bssn 
obserted  in  wind  tunnel  sxperiiwnts  and  free  fli^t  tests.  On  the 
other  hand,  conclusions  drawn  fron  stability  thsory  were  eonpistely 
contrary  to  obsenratlon  on  blunt  objects.  A  considerable  farorahle 
pressure  gradient  usually  exists  on  such  objects  orer  the  region 
which  is  coTsred  by  boundaiy  layer  flow,  and  fron  this  fact  one  expects 
the  transition  to  be  delved,  fixperinents,  bowerer,  showed  that  tha 
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transition  rieynclds  nunber  based  on  distance  x  is  of  the  sane  order 
of  Maeni  tjidtt  as  03  s  flat  plate  under  Isotbencal  flow  coaditions. 

Cool!  jg  of  the  surface  was  found  to  have  very  little  effect  on  the 
transition  Btynolds  nuaber.  Seference  367>  for  instance,  reports  tbs 
results  of  experinents  on  sphere-cylinders  and  cn  elliptic  cylinders 
at  a  ratio  of  total  enthalpy  to  wall  enthalpy  between  9  and  30  and 
at  a  Mach  saber  around  3.  The  experinents  were  perfonad  in  a  shodc 
tabs,  and  it  was  found  that  the  transition  oeearred  at  a  H^molds 
saber  based  on  nonentua  thickness  between  225  and  325  (the  last 
fignre  for  an  enthalny  ratio  equal  30).  Ibe  transition  iiqrnolds 
sster  based  on  length  x  was  between  500,000  atxi  10^.  Fires  the  above 
figures  it  can  be  seen  that  cooling  had  no  pronounced  effect.  Qa 
shape  InfluenctKl  the  transition  itoynolds  ntaber  sonewhat. 

another  observation  for  which  a  definite  explanation  is  still 
sissing  is  the  fact  that  on'  slander  objects  tuxbnlent  flow  has  bean 
obeartad  at  extresely  strong  cooling  of  tha  surface  and  under 
conditions  where  the  stability  thaory  predicts  conpletely  laninar  flow. 

Ihs  nonsotun  thickness  as  a  funeti^  of  the  distance  x  has  to  be 
known  if  ilsynolds  nuabers  based  on  nossntun  thiclaess  are  to  be  used 
for  a  prediction  of  the  location  where  the  laninar  boundary  layer 
ehangee  into  a  turbulent  one.  Calculation  of  the  aosantus  thickneas 
is  considerably  sore  tedious  than  that  of  heat  transfar.  In&efersnee 
43!),  Lees  proposes  an  iteration  procedure  by  which  the  mcmimtm 
thickness  can  be  calculated  for  objects  of  arbitrary  shi^s. 
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/yg.  X? 

Vlscosl'.y,  ent  ccnduotivity  an'!  soecV;  *  hed»  o:'  -.ir  at 
1  ata.  pressure  «5d  in  dissociated  cqu  -ibriua  'from  Ref,  U27) 
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7.  t»aoFamES 


A  knovledgs  of  the  thsimodynsMio  end  transport  properties  is 
required  for  the  gases  participating  in  the  flow  when  one  wants  to 
evaluate  an/  of  the  relations  contained  in  this  report,  for  air,  these 
properties  have  been  calculated  in  a  large  pressure  range  and  up  to 
vex7  high  tenperatures  using  the  methods  of  statistical  aeebanies. 

Such  infbxmation  is,  for  instance,  contained  in  Beferenees  ItH  and  1(27. 
Oal/  indirect  verification  is  available  b/  the  fact  that  beat  transfer 
coefficients  ■assured  in  shock  tubes  at  very  high  teaperatures  agree 
quite  well  (within  13-20  %}  with  predictions  based  on  the  relation* 
discussed  in  the  previous  sections  and  on  the  calculated  air  properties. 
Tisoosit/,  heat  oonduetlvlt/,  and  specific  heat  of  air  under  conditions 
close  to  tbezvod/nsBic  equilibriun  are  presented  in  Pig.  13  for  one 
atsu^pbere  pressure.  The  values  have  been  taken  fromaffwretwj;  ii27» 

Pig.  lii  contains  the  itl'sndtl  mnber  and  the  lewis  zsafoer  desoriUng 
diflhsion  of  dissociated  aUws  again  for  air  at  one  atnosphere 
pressure,  it  has  been  pointed  out  before  that  two  definitions  are 
in  use  for  the  s-^ifie  beat  and  for  the  beat  conduetivitgr.  According 
to  one  definition,  the  specific  heat  based  on  the  enthalpQr 

which  contains  not  ocl/  sensibl*  heat  but  also  the  chaaioal  energ/. 

The  second  definition  bases  spseifie  bsat  (c^)  on  tbs  ssnsibls  bsat 
oni/.  In  s  sinllsr  wa/,  ths  bsat  eonduetlvit/  can  be  defined  b/ 
oonaldering  the  transport  hgr  heat  conduction  alone  (k) ,  or  bgr 
including  the  energ/  transpert  caused  hr  diffusion  (k*).  Ths  second 
part  is  uniqusl/  related  to  the  tenperature  and  pressure  since  it  is 
postulated  that  conditions  near  equilibriun  exist  in  the  gee.  Pig.  13 
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Fig.M 

Prandtl  ".uaber  an'  Lewis  number  ','f  air  at'l  atm 
and  in  dissociated  equilibrium  (from tel',  lt27^ 
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shovs  that  the  difference  between  the  Talues  of  specific  heat  and  of 
heat  ccnductiYity  according  to  the  two  definitions  is  very  large. 

The  values  according  to  the  definition  which  considers  sensible  heat 
only  have  to  be  used  in  the  relations  presented  in  this  report,  ilie 
tVandtl  and  lewis  nuobers  in  i'^Lg.  Hi  are  based  on  the  heat  uunductivitor 
k  and  on  the  specific  heat  Cp. 

In  addition,  the  connection  between  enthalpy  and  tei^>erature  has 
to  le  known  for  calculations  of  heat  transfer.  For  a  gas  in  thermal 
equilibrium,  this  connection  is  best  obtained  from  enthalpy  entropy 
diagrams.  (hoUier  Diagram  for  Equilibrium  Air,  Aveo  Assearch  Lab., 

Avco  hanufacturing  Corp*,  Sverett,  Hass.,  Jan  1957i  HoUier  Chart  for 
Air  in  Dissociated  Equilibrium,  KAVOdD  riep.  UltU6,  liaval  Ordnance 
Laboratory,  Silver  brings,  Hd.,  1957,  see  also  Ref.  Ii39). 

The  situation  is  considerably  more  involved  to  cp'  rul.-.t-'ons  on 
mass  transfer  cooling  processes.  The  detemlnation  of  the  properties 
involvea  in  this  process  is  usually  the  most  time-consuming  prerequisite 
for  such  calculations,  especially  when  chemical  reactions  occur 
within  the  boundary  layer  as  they  have  been  discussed  in  Section  L2. 
Information  on  the  properties  of  the  vaurlous  constituents  involved  in 
mass  transfer  cooling  is  usually  quite  restricted,  and  the  calculation 
of  the  properties  of  the  various  mixturos  which  is  possible  by  the 
methods  ofstatlstical  mechanics  is  very  involved.  The  reader  has  to 
be  referred  to  the  literature  for  the  detemlnation  of  such  properties. 
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Ai>Pia>Dii  I 


Use  of  teaper&cures  or  entti&lpies  in  best  transfer  calculations 
has  an  inO.nanee  on  the  energy  balance  dascribing  heat  exchange 
betveen  a  wall  and  a  nere  component  sixture.  This  will  be  denonstrated 
in  the  following  exanple)  which  considers  flow  over  a  mass  transfer 
cooled  wall,  llte  coolant  nay  be  denoted  by  1.  and  it  will  be  assunad 
that  the  wall  surface  is  inpemaable  to  the  gas  2  in  the  nain  flow, 
let  us  now  fix  attention  to  tbs  mass  flow  of  both  coaponents  1  and  2 
through  a  plane  a-a  arranged  in  the  gas  Just  outside  and  parallel  to 
the  surface,  llw  ness  release  gMsrates  a  conweetise  relooi^  t 
through  this  plane  and  «.  mass  flux  of  coaponent  1  equal  lAere 

6 .  is  the  partial  density  of  1.  Ihere  is  also  a  nass  transport  by 

in  . 

diffusion  of  snount  -  ^0(  'w*  ^  indicatss  tha  nass  fraction 

w,=  ii.  Tha  total  transport  of  conponsot  1  is: 

y  dw. 

f  V  -  ^w  *  “ 


fbr  conponent  2,  tha  nass  flow  throng  the  surfacs  and  tharefore  also 
through  ths  control  plana  a«a  ia  xaro: 

dw. 


fV-  =  o 


(1-2) 


Mow  wa  consider  ths  energy  transport  through  ths  control  plans,  lha 
energy  flux  per  unit  tine  and  area  nay  ha  denoted  by  q* .  According 
to  equation  (3-38) >  the  energy  flux  ia: 

q'  =  -k(  ■  fWi(  4?^  ‘ 
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This  energy  flux  must  be  equal  to  the  flux  through  control 
plane  b-b  arranged  in  the  solid  wall  ;}ast  beyond  the  interphaae.  The 
Ints'-vhase  moves  continuously  towards  the  intr;ricr  of  the  wall  because 
of  the  mass  release,  k  mass  flow  of  solid  waterial  with  the  ness 
Telocity  n  exists,  therefore,  through  the  plane  b-b,  which  is  considered 
fixed  relatiwe  to  the  interphase.  A  corresponding  energy  flnx  by 
convection  is  nig,  where  i^  indicates  the  enthalpy  of  the  solid  aaterial. 
Another  energy  tr.»nsport  by  conduction  nay  be  indicated  by  q^. 
Conservation  of  energy  requires: 

q*  =  qg-*'  fflig  (I-U 

Introduction  of  equation  (I-l)  into  (1-3)  results  in: 

S,  =  -  -37  -  V 

Vie  have  here  obtained  the  result  which  was  expressed  in  c’  aptar  2  by 
equations  (2-3)  and  (2-9)  with  the  exception  that  no  radiation  was 
considered  in  the  present  derivation. 

The  heat  flux  q,  aay  now  be  expressed  by  the  enthalpy  gradient  in 
the  gas  mixture.  With  n~^r  end  equations  (1-3)  and  (I-U)  one  obtains: 


The  first  tern  in  this  equation  is  idr«ntical  to  equation  (3-36),  since  / 

^  Owi  dw2  ; 

in  a  two- coBpcnent-aixturo  -gy  =  "  and  since  -  i2*  I 

For  a  gas  mixture  with  1«  =  1,  the  first  ten*  takes  the  shape  of 

ecuatlcn  (3-4:2;,  and  the  heat  flux  q^  is: 

Si'- --r  V 

P 

i  is  the  total  enthalpy  of  the  gas  mixture,  1^  indicates  it  value  at 
the  Interphase,  and  i^  is  the  enthrlpy  of  the  ablating  'satorlal  at.  wall 
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surface  tenperature  and  In  iiolld  ^tate. 

in  equation  on  the  other  band,  the  enthalpy  1^  of  the 

ablating  material  in  toe  gaseous  stats  and  at  wall  surface  temperature 
takes  the  place  of  i^  in  equation  (1-6).  The  enthalpy  i^  can  be  more 
easily  determined  in  most  applications  than  i^,  and  the  form  (1-5)  of 
the  beat  flow  equation  is  therefore  preferable.  In  the  use  of  beat 
transfer  coefficients  for  mass  transfer  cooling,  one  has  co  watch 
carefully  whether  the  coefficient  was  defined  by  equation  (2-3)  or 
by  a  defUiition  .:hich  replaces  the  temperature  gradient  in  this  equation 
by  an  enthalpy  gradient  as  aigfested  by  equation  (1-6).  For  a  gas  with 
properties  depending  on  temperature  only  there  is  no  difference 
betMeen  i^  and  i^,  and  both  equations  (1-5)  and  (1-6)  becoms  identical. 


7 


‘O' 
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The  mass  iraction  of  aiy  of  the  coapooents  in  a  tHO~coaponent 
ailxturr.  J.S  determined  a  diffusion  eqoation: 

Oyi.  d  wi  0  d  lb 

■  ■5y  “i 

denotes  the  production  of  the  species  i  bgr  the  chemical  raaci  ica 
per  unit  tine  and  Toluiee4 

The  diffusion  equations  can  also  be  written  in  ‘•psewlo-aess- fractions" 
w^.  A  pseudo-aass'fraetion  givcs  the  mass  fraction  of  an  element  in 
the  aixuire  regardless  in  what  compound  it  is  contained.  The  pseudo¬ 
mass-fraction  of  oxygen,  for  instaaee>  is  the  ntio  of  the  orqrgen  mass 
to  the  total  mass  of  the  mixture  regardless  whether  osgrgen  occurs  as 
0,  OO2  or  in  aigr  other  compound,  so  element  is  destroyed  or 
created  in  a  chemical  reaction.  Tharefore,  the  diffusion  equation  in 
pseudo-mass-fractions  reads: 

dvT  3?:  0  , 

+.  y  -^)  5  ^  )  (11-2) 

Fbr  a  gas  with  Pr  =  1  and  le  s  1,  equation  (11-2)  has  the  same 
form  as  the  energy  equation  (i(-8).  This  fact  can  be  utilised  to 
calculate  the  composition  of  the  gas  throughout  the  bouDdaxT’  layer  and 
a'  the  wall  surface.  L.  lees  carried  this  idea  through  in  Tarlous 
exanoles  (nef.  1*32). 
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